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CHAPTER - I 


GEHERAL IHTROHJCTION 

1,1 IITTROPJCTIOW ; 

Thou^ the iise of lubricant in relatively moving parts to 
facilitate motion is known to man kind since the invention of yiheel 
and axle, it was only in 1886 that the real mechanian of lubrication 
was explained by Reynolds [1886 ] in his classical paper. The 
obj ect of lubrication in mechanical system is to minimize wear 
and enerQT loses by reducing friction between the two relatively 
moving surfaces (plane/curved) by the presence of an interposed 
film of lubricant, 5he proper knowledge and understanding of the 
lubrication process are essential in in^jroving the standard of design 
and the efficiency of the system. The characteristics of the bearing 
system depend upon ihe following types of lubrication process, 

1, BEydrodynamic lubrication 

2, Boundary lubrication 

3, Mixed lubrication. 

In l^drodynamic lubrication, the sliding surfaces are separated 

-5 -3 

by a layer of fluid -film of thickness of the order 10 - 10 cm and 

the load capacily and the friction farcre depend upon the viscosily 
of the lubricant and the gecmetry of the film, Pinkus and Stemlicht 
[1961 ], 
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In bcHindaiy lubrication, film-thickness between liie two 
surfaces is veiy small (of the order of a molecular thickness) 
and Uie friction is determined by the physical and, chemical properties 
of liie film and the solids in contact. Boundary lubrication is 
assoaiated with the effect of adsorbed iiydrocarbon molecules vdiich 
seperate the two surfaces and the real contact takes place along 
the *tails* of molecxxles sticking out of the surfaces j Bowden and 
Tabor [ 1950, 1964 ] and ling et. al. [ 1969 ] . 

In the case of mixed lubrication, a situation between hydro- 
dynamic and' boxindaiy lubrications, the fluid film between the two 
sliding surfaces is thin enough and the surface asperities begin 
to interfere in the hydrodynamic process. Here the load is partly 
supported by asperity contacts and partly by fluid film.' Ihe force 
of friction also depends upon the interaction of the surface 
asperities and the action of the film, lUller [ 1954 ] , 

Lenning [i960 ], Dobry [ 1964 Christensen [ 1972 ] and Tsao and 
Tong [ 1975]'. 

In gener^, the lubrication characteristics of the beai^g 
system depend upon the followirg; 

(l) Hature of lubricant t liquid or gas, its Hewtonian or non- 
Hewtonian behaviour, variation of viscosity or consistency with 
temperature, pressure or concentration ^diidi may change along a*; 
well as across the film-thickness, its affinity with surface such 
as adsorption, absorption, chemical reaction etc,, Braithwaite C1967J:; 
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(2) Uature of surface; surface roughness, elastieityt thermal 
conductivity, hardness, porosity etc., its affinity with the 
lubricant. 

(3) Effects of flow regimes , bearing configuration, film-^thickness, 
boundary conditions and the surrounding environment etc. 

!Eo study the behaviour of any lubricated ^stem a mathematical 
model is made by taking the above factors into accoiint whidi depend 
upon a given physical situation. dhe bearing characteristics such 
as load, dflow flux, friction force etc, depend upon the pressure 
generated in the film, and the lubrication process. Ihe surv^ of 
varfous attempts, made by many werkers to obtain the governing 
equation for pressure in the lubricant film, are summarized hei^e, 

1.2 sCTEVEf OF REgNoiDs aagAfioij 

In the case of hydrodynamic lubrication, the equation governing 
the fluid film pressure in a bearing can be obtained by coupling the 
equations of motion and the equation of continuily, lEhis equation 
was first derived by Reynolds [1886 ] in his classical paper by 
considering following assiaaptions; 

(i) Ohe radious of rairvSi'i^e of the bearing surface is 
large in comparison to the film-thickness. 

(ii) Bae lubricant is an incompressible Newtonian fluid, 

(iii) Bie lubricant viscosity is constant- 
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(iv) Inertia and. body force terms are small incongjarison 
to the viscous and pressure terms » 

(v) Owing to the anallness of the fluid film-thickness, 
the velocily gradients across the film are large 
incomparison to the velocily gradients along "ttie film. 

(vi) illiere is no slip at the fluid and solid interface. 

Since early fourties, various attoapts have been made to 
generalize the Reynolds equation applicable to bearing ^stem 
functioning under unusual conditions such as hi^ temperature, hi^ 
pressure etc. 

33ae fifst attempt in this direction was made by ihgg [ 1946 } 
w4io proposed the themal-wedge concept in the lubricant film. 

Cope [1949 1 relaxed the assumption made by Reynolds [ 1886] to 
extend the theory for variable viscosily and densily along tiie film. 
But the variations of fluid pressure and fluid properties across the 
fluid film were still neglected. In [1950 ] , Wannier has shown 
that the basic Reynolds equation can be derived from the eouations 
of motion by considering the variation of fluid pressure across the 
film. Later, Halton [ 1958 3 extended the investigation performed 
by Cope [ 1949! I and derived a generalized form of R^nolds equation 
f<m* rough bearings. Several other studies lave also been conducted 
by considering the temperature and viscosity variations along as well 
as across the film, Zienkiewicz [1957] , Cameron [ 195& J, , Cameron and 
Wood [ 1958] fHuntex^aad. ‘.Zienkiewicz [196(1], and Yoxang [1962 ] , 



5 


In early sixties, an unified approach was made by Bowson [1962]' 
to generalize the Reynolds equation by considering the variation of 
fluid properties across as well as along the lubricant film, Dowson 
and Hudson [l963-a, l963-bl f Curie et. al, [1965] t Hahn and 
Zettleborou^ [1967] j McCallion et, al. [1970] have also studied 
the effect of viscosity variation in lubricated systems under isothermal 
or adiabatic conditions by using numerical techniques. Recently, 
the performance of viscosity variation across as well as along the 
hydro dynamically lubricated film was studied by Quale and TIH.ltahire 
[1972] , Ezzat aM Rohde [ 1973 -a, 1973-bl . 

A different approach to study the effect of viscosity variation 
has also been proposed by ELpei [ 1962 ] , iCLpei and Hica [ 1967] and 
(ELpei and Degueurce [1974] by assuming a relationship between viscosity 
and fLlm-thiclcness for convergent films, 

SUP EPFECfS: 

In the above mentioned studies, one of the b^ie assumptions 
used is no-slip condition at the bearing surface. However, in 
general this assumption is not valid and the effect of slip may be 
important on the flow behaviours of gases and liquids especially 
when the fluid film-thickness is very small, Kennard [ 1938 ] , Schaaf 
and Chambre [ 195& ] } the surface is very siaooth, Devenport [1973] ? 
and when the bearing surface is porous, Beavers and Josejh [1967] , 

In the case of gas lubricated bearings, as the density is lowered. 
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■fche molecular mean fi?ee patit becomes comparabj.^^^;..^ *411® ■tdaicknessj 

■fcbe gas seems "to lose its grip upon tbe surface and leaves it with 
a finite slip-velocity, Kennari [1938^ , Srad [1949] ani Schaaf and 
C9iambre [195S ] . 

Birgdorfer [ 1959 ] was the first to derive a modified form of 
Heynolds equation applicable to hydrodynemicalHy lubricated gas 
bearings with slip— flow condition, Gross [ 1962 } . Hsing a nd 
Malanoski [ 1969 ] have studied the effect of molecular mean free path 

in the spiral-grooved thrust bearing by considering the ^ip-velocity 

* 

at the surface* !Ihe rarefaction effects of gas lubricated bearings 
in magnetic recording disk file has also been studied, Tseng [1975] . 

In the case of liquids, the slip flow can be in^Kirtant when the 
surface is smooth and bearing is functioning at high temperature, 

Iamb [1945] , Bird [1958] and Bevenport [1973] . Slippage at the 
wall is more in the case cf high temperature where the ■viscosity of 
"the base oil decreases near the surface* This effects may be eliminated 
by roughening the surface, Bramhall and Hitton [i960] and Bevenport 
[1973] . 

The existance of slip flow at the porous boundary also been 
demonstrated by Beavers et. al [1967, 1970] and. it has further been 
supported by Saffiaan [1971] and Taylor [1971 ] . The effect of slip 
flow in porous bearings has beoa studied iy Sparrow et.al [1972] , 

Wi [1972] , Jones [1973] , Murti [1973] and Beavers et, al.[ 1974] • 
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In general, the effect of slip is to decrease the load capacity 
of the hearing and the friction force at the hearing surface. 

ROUGHHESS EPPECTS 

As it has been pointed out earlier, the topology of the hearing 
surface plays an important role in the lubrication mechanism. In one 
of the approaches, the effect of surface roughness is taken into accoun-t 
in the usual Reynolds equation by considering that the film-thickness 
is a function of surface roughness which may he represented by a series 
of sine or cosine waves, Burton [1963I , Shukla and Prasad [1966] . 

This procedure has also been applied to study the bearing characteristic 
of rollers, Dowson and Whomes [1971] and of spiral grooveaheating, 
Wildmann [1968] . 

In another method, the surface rou^ness is assumed to be 
represented by a stochastic process and the usual procedure is followed 
by taking statistical mean or average of the basic governing equation, 
Papoulis [1965) • Using this approach, the characteristics of an 
infinite slider bearing and short journal bearings have been studied 
by considering the distribution of the surface roughness to be gaussian, 
Tseng and Saibel [l967-a» 1967 -b] . Using this approach, Christensen 
[1969-70 ] has derived the generalized form of Reynolds equations appli- 
cable to finite rough bearings and its various particular oases have 
been investigated, Christensen [1971} , Christensen and Tender [1969*1971 
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louder aud Cliristensen [l972-«,, 1972-b] . Arioliier form of Reynolds 
equation S 5 )plieal)le to rou^ or deformed surfaces was presented by 

: , ' ll 

Berthe and Godet [ 19731 • A further refinonent of this Idaeory has 
been made recently Christensen et^al, [1975] • The concept of 
stochastic approach has also been extended to rou^ beairings using 
compressible lubricants, ELord [1973] * Christensen, et* al* [1975] . i 

1.3 MIXED HJBRICATIOH' 

In the above, a airv^ of Reynolds equation applicable hydrodyna- 
mically lubricated rough beaiii^gs where the asperities heights are 
very small incomparison to the rKjrminal film thickness has be^ 
pi^ented. However, ?iien the heights of the surface asperities are 
of the sane order of magnitude as the fLlm-thickness the direct 

■ I'. 

contact between the opposite surfaces may occur, causing liie coMition 
of mixed lubrication. In this case, the total load may be partly 
supported ty the asperities contact and partly by hydrodynamic action, 

I 

Similarly, the total friction force will arrise partly from boundary | 
lubrication condition associated with the asperities contact and j 

partly from tydicdynamic action, lUller [ 1954] , lenning [i960] 
and Dobry [l964] . C9aristensaa [1969] and lowles [1971] have 
studied the asperities deformation and hydrodynamic effects in 

i 

lubricated contact, Bie transition from boundary to mixed lubrication 
has also been investigated by Hiomson and Ulliam [1972] , Recently, 
a -yieoiy of mixed lubrication by usir^ stochastic R^nolds equation 
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has been presented by Christensen [ 1972) . !I?sao and !Ibng [1975] 
have also shidied the condition of mixed lubrication by calculating 
1he roxi^iness interference ecrea frcm •te knowledge -of bearing 
geometoiy and surface roughness. 

Keeping the above in view, in this thesis an attoapt has 
been made to give new ideas and to generalize some of the older 
theories preserved before in the theory of lubrication* Ihe 
following summary gives the details of the work done presented in 
•fcis thesis. 

1.4 SDMMAIg 

Dais thesis consists of five chapters. Chapter-I deals with 
the general introduction of the literature related to 1he work 
presented in the thesis. 

In Chapter-II, a generalized form of Reynolds equation wilh 
slip at the bearing surfaces is derived by considering the variation 
of flviid properties across -as well as along the film. Itor gas- 
lubricated bearings ,a generalized form of H^nolds eqviation is 
deduced in terms of molecular mean free path. Hie case of hydrostatic 
bearing is studied and it has been shown that the load capacily 
decreases as moleculajr mean fire^ path increases. 

In the case of liquids, a generalized form of Reynolds equation 
is derived by considering the viscosily of the lubricant as a step 
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ftuajtioii across the film in nnilttple l^er lubrication. aiie 
case of slider and hydrostatic bearings are studied. It has 
been shown that the load capacity^ friction fbroe increase as the 
coefficient of sliding friction and viscosities of layers increase. 

lUrther, in the csise of porous bearing, a generalized form 
of R^nolds equation is deduced by assuming the existance of slip 
flow at the porous surface. Moreover, the c se of one dimensional 
externally p3?essurized porous squeeze-film bearing has been 
investigated. 

In Chapter-Ill, by using the theoiy of stochastic process, 
a generalized form of Reynolds equation applicable to rou^ bearing 
is derived by assuming that the fluxes axe represented by poww: series 
of stochastic fllm-thidmess function, OSie case of hydrodynamic 
one^imensional step bearing and hydrostatic bearing are investigated. 
In the case of a step bearing, it has been shown that the load 
oapacily and friction force increase but the coefficient of friction 
decreases as the rou^tiness parameter increases. Similar results 
are true in the case of a hydrostatic bearing, 

CSbapter-I? deals with the generalization of the vroife presented 
in Ghapter-III, ty (considering the variation of viscosily alor^ as 
well as across the fLlm-thickness and stochastic form, of Reynolds 
equation >has been derived. The effects of viscosity variation in 
■file cases of slider and tydrostatic bearings have been studied. 
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In Chap ter -Y, a new theoiy of missed lubrication is presented 
try eonsideriog the surface asperities to be rand<mily dritribw<t«d fins 

or cxrnes idiich peneratrato thrcHi^ Ihe fluid film and j^rilly may be 
in contact with the asperities of the opposite aiif ace. !Ili^ are 

considered to foim a net-woife of interacting aone corresponding to 
each surface throu^ which the lubricant flows as happens in ihe 
case of flow through porous channels with singular resistances or 
some Tidiat as in the case of flow throu^ porous matriac* 

A generalized form of Reynolds equation is derived by considering 
a modified foms of ITavier-Stoke’s equations in Ihese interacting 
zones. As pariicular cases, step slider bearing and hydrostatic 
step bearing have been studied. It has been shown Uiat the load 
capacity and friction force increase as the interference to 1he flow 
increases or as viscosities of layers increase. 

Ihis theo3y has also been ^plied to ^e study of synovial 
joint. It has been pointed out that the load capacity and time 
of squeezing decrease in the case of fractured joint under impact 
loading condition^ ifowever, this my be c<sapensated by boosting 
effects which increases the viscosity of the i^novial fluid. 
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CHAPTER - 2 


EFFECTS OE SLIP AKD VISCOSI5Y VAEIAHON IH HTDROOTAMIC IDBRICATIOl 

2.1 IHTBCIOQTIOM 

In general, most of the lubricated systems can be considered 
to consist of moving (stationary) surfaces (plane/curve, loaded/ 
unloaded) with a thin film of an external material (lubricant) between 
them. Hie presence of such a thin film between these surfaces not only 
helps to support considerable load but also minimizes friction. Hie 
characteristics, such as pressure in the film, frictional force at 
the surface^ flow rate of lubricant etc. of the system depend 
the nature of the surfaces, the nature of the lubricant film, boundary 
conditions etc. 

Hie equation governing the pressure generated in the lubricant 
film can be obtained by coupling the equations of motion n ryi the equation 
of continuity. Hiis equation was first derived by Reynolds [1866 } in 
his classical paper under usual lubrication assumptions for an incom- 
pressible lubricant and is now known as 'Reynolds equation'. In 
deriving this equation, the thermal and compressibility effects as 
well as viscosily variation in the film were ignored. Later Cope [1949] 
modified the R^nolds equation by including viscosity and density 
variation along the fluid film. Hie viscosity variation across the 
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film thickness has been, considered by Zienkiewicz [1951,1957 ] , 

Cameron [1958 ] who also pointed out that temperature gradients and 
viscosily variation across the film may not be ignored. In the year 
[1962 ] f Dowson unified the various attempts in generalizing the 
Reynolds equation by considering the variation of fluid proiperties 

i 

across as well as along the fluid film-thickness by joe^ecting qhip | 

effects at the bearing surfaces. 

nae effects of viscosity variation have also been considered 
by Qttaie and Wiltshire [1972 ] and lipei and Degueurce [1974] without 
considering slip at the surfaces. However, it has been noted that the 
effect of slip may be injiortant on the flow behaviour of gases and 
liquids particularly when the film-thlekness is very Knall, Kennard 
[1938 ],Srad [1949] and Gross [l962 land the surface is very smooth, 
Devenport [1973 ] . fhe effect of slip can also be important at the 
porous boundary, Beavers and Joseph [1967] . 

In the gas bearing applications, such as gyroscope, Heinrich 
[l970] , the bearings are essentially operated in the slip-flow regime. 
Ihis happens because of either low pressure or extreiiely thin bearing 
film when the length of the molecular mean free path X becomes comparable 
to the film-thickn^s h, Kennard [l938] and Grad [l949 ].» Bor smaller, 
but non-negligible value of ^h, the hmediate adjacent layer of the 
gas to the solid , surface is no loiter attached to the surface but has 
a finite relative slip velocity and this produces the effect of an 
apparent dlminuation in the viscosity of the gas. 
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3b study the effect of slip, Burgdorfer [1959] modified 
Eeynolds equation for gas-lubricated bydrodynaBiic bearings mth 
’’slip flow” under isothermal condition. He pointed out if 

0 ^ < 1 , gas flow may be considered to be continuous and analysis 

can be carried out with modified slip boundary conditions, Also 
Hsing and Malanoski [ 19691 studi^ the effect of molecular mean free 
path in the spiral-grooved thrust bearing and 3bing [1975J used the 
R^nolds equation to study the rarefaction effects of gas-lubricated 
bearing in magnetic recording disk file. 


Slip flow boundary condition at the wall in the case of 
gas bearing can be written as Kennard [ 1938] and Burgdorfor [l959] , 


u 


slip 


* 0 ( 5 ^) 


X(^) + 

^az'^wall 


(2.l) 


where f is reflection coefficient, X is the mean free path and o 
is numerical constant. Because of_a and f are closed to 
unily, it mgy be considered that ^ (""J") is unity. As Molecular 
mean free path X depends upon the fluid viscosity, pressure and 
temperature it can be approximated by the following relation, 
Tsing[l975]: , 

Xe-^ - /rt (2.2) 

where R is a gas constant and 3? is the temperature, n is the 
viscosity of the gas and p is its pressure. 
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She effect of wall slippage in the case of liquids witti 
dispersed mateidal and greases have been discussed by Devenport 
[ 1973 ] * Slippage may be more at higher temperalwre \idien the 
base oil viscosity reduces but this m^ be eliminated by rou^ening 
the solid boundaiy. Hie magnitude of slip depends upon the type 
of fluid and surface roughness. Mathematically, the effect of 
slip has been pointed out by Lamb [ 1945 1 , Bird [ 1956 ] and 
Bramhall and liitton [ i960 } , Hie velocily of slip at the wall in 
the case of liquid, can be defined as follows: 


%lip 
inhere 6 


_ JL / iii') 

“ e 9z^wall 

is -the coefficient of sliding friction at the wall. 


(2.3) 


Hie slip flow is not only iii^orfeant at the solid bearing 
surface but it plays an important role in bearings with porous 
facing, Hiis effect has been discussed by Beavers, et. al. 

[1967, 1970 3 for an incompressible fluid in their eicperiment vfco 
demonstrated the existance of slip-veloci1y at the porous surface 
and proposed an aaperloal boundary condition. Hie work was further 
supported by Hay lor [1971 ], Saffinan [l97l3 by providing theoretical 
justifications. 


Effect of slip flow in liie case of porous squeze film and 
porous spherical shell were considered by Sparrow, et. al. [1972] , 
Wu [ 1972 ] , Jones [ 1973] , and Marti [1973] . Beavers et. al. 
[ 1974 ] have also proposed a boundaiy condition at a porous surface 
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and experiments were perfomed to establish the slip flow for porous 
gas bearing. Glie slip velocily at the porous surface can be 
written as, Beavers et. al. [1974] 



(JH) 

'■82'^ wall 


(2.4) 


where § is the slip coefficient at the wall and 4 is the pemeability 
of the porous facii^. 


In view, of the above, in Ihis chapter, a generalized fcrm 
of Reynolds equation for fluid film lubrication wilh slip velocity 
at the surface is derived, by considering the variation of fluid 
properties across as well as along the film-thickness, and various 
particular cases are studied. 


2,2 BASIC IQUATIOITS ; 

, I 

Consider the laminar flow of a fluid betweeh tsro 
surfaces T«hose pt^rsical configuration is shown in rig, (2.l). 
Considering the variation of fluid properties acsross as well as 
along the film-thickness, the basic equation of motion and equation 
of continuity in their most general form for a Newtonian fluid can 
be written as follows: 




ax 


.i2) +1 1. . JE) 

ax ay 5 ax i ax 
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pY.J£ + 2 Jl. 1- n(^ 

ay 3 ay ay ax'^ 3 ay ^ay 


, ° . av\ . d /<?V , dU\ 

3z ’^^ay az^ ax ’^^ax^ ay^ 


' av , au ' 


(2,5 ) 


p7_ii£ + 2^ +-^-^ nfiS-il) 

^ az 3 az ^Sz ax^ 3 az ''az ay^ 


3 nris 4. J. i- 

ax ^az ax ay ay az 


ap ^ j_(pu-) i_(pv) ^ j__(pw) _ 
3t ax ay az ~ 


( 2 . 6 ) 


To derive Uae generalized form of Reynolds eqviation vd.th 
the usual assumptions of lubrication theory are made,I)owson [l962] 
following him, the equations (2.5) could be simplified as follow: 


jR = 4_ (n -M) 
ax az '■ az'^ 


(2.7) 


ap = 1- (n il) 

^ az ^ az^ 

where p = p(x,y) is the pressure in the film. 

Considering the slip at the surfaces, Ihe boundaoy conditions 
can be written as follows: 


= + V. , z = H ’ 

Vaz^i 1 ^ 1 


( 2 . 6 ) 


U = (u)^ = - 'Tj- “ = 
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where X's and 6's are molecular mean free path for gas lubrication, 
and depend upon lubricant temperature, pressure and viscosity. In the 
case of liquids they depend upon viscosity and the coefficient of 
sliding friction- However, in the case of porous bearing, these are 
functions of slip coefficient at the v/all and peimeability parameter 
of porous facing, [see equations (2.1) to (2.4)] . 


Integiating equation (2.7) and using boundaiy conditions 
(2.8), the following expressions for the fluid film velocities are 
obtained. 

z 3P ) I'. BP 2 , 

u = IJ, + [a.H, J — ] — + [ ^ ^ ] K + i “ ] 

1 114 n ax p pax i 4 n 


where 


H. 




ap 


dz 


H 


2 dz 




H. 


(2.9) 


H, 


P. 


1 = “l"*! * V2 * /. 


2 z dz 
H, n 


(2.10) 


H. 


PJ 


= 8 H + 6 H + r 

1 ^1 1 ^ 2 ^ I 


“1 n, ’ “2 >2 

*1 .^2 
1 ° '>1 ’ ■ "2 ■ 


2 zdz 
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Integrating equation of oontimilgr (2.6) with respect to z 
and taking limits fcccm z = to z = we get 

jj jj 

^ (p^) dz + / ^ — (pv) dz + lpw3g^ - 0 




Applying the well known integral formula 


(2.11) 


A 


H, 


^ I ^ f(x,y,z) dz = /' ~ {f(x,y,z)} dz + f(x,y,H ) 


3H„ 


H. 


H, 


2' 3x 




-f(x,y,H^)-^ 

to the equation (2.II) and siB5)lifying, we have 


(2.-V2) 


H 


^ h ^ b - <->3 ^ 


H. 


% 


ay 




9H 3H H 

+ (Pu)^ ^ + (pt)^ ^ + [Pw]g2 ^ Q 


(2.13) 


The integrals of (Pu) and (pv) can be evaluated by integration by 
part to give the following equation : 

I 1% dz + tlj (pu)^ + (pv)^] -H, [-^(pu), ^--|■(pv),] 


fet I ^ If* ~b f / ^ 


H. 


H. 


9z 3z' 


H 


* fpw3g^ = 0 


(2.14) 



Introducing the expressions for u,v and their derivatives in 


the above equation (2. 14), we get, 


3P 


3P 




= ®2 (“"'a * W 


8y 


[ jT +n, 




- (Vp - ^i) (^'•5 + (>p) 

l~ +7, G, ] 


3y 


^o 


’1 ^3 


H 




H 


H. 


wheare 


J'P = / 


^0 .. 

— f ^2 


H. 


11 

P 


) dz 


H 


Pi 


H. o 


S’, = J 


H. 


2 p z 


H. 


dz 


S, - t z -r- + j' — ) - ^ 


‘ H. ^ 0 ' H. ^ 


1 


o 

- p» 

1 

u 

0 


1 


(2.15) 
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H 




dz. 


H. 


1 


Equation (2.15) represents a generalized form of Reynolds equation 
-for a fluid film-lubrication wlaeTe slip-velocities are considered at 
the bearing surfaces, The two sets of functions P and ff depend upon 
the variation of fluid properties along as well as across the film 
and slip conditions at the surfaces. 

In the case of no-slip at the boundaries i. e. 


S ='‘2 = ° 

“l =“2=®1 -82 = ° - 

the generalized equation (2.15) and equations (2. 16) reduce to the 
following forms 

9P 




= (p u) + ~ (p v) ] -H. [■~(ipu) +|— (p v) ] 

^ 3x 2 3y 2-' 1 3x 1 3y 'I-* 

„ ~ - A .. . . .. . 4- U cl] _ -S— r , . . 4 .— J— . . . -i? — . M. 4. y 0- f 

ax'- ¥ I ^ 


3P , ^2 

®1 ^ 


(2.17) 
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vshere 


H 

jt _ f 2 ^ 

O ^ yk 




f ^ z <iz 
^1 ' 


H. 


S'. 




H. 


1 


(2.18) 


H 


P_ = / 


H. 


2 P 2 dz 
n 


1 


=/"^ 

", 

®2=J 

", 


<^ - / 

T1 








H. 


r (r 

i; V 


H 


1 


z dz 


T? J 


O K 


dz 


)] dz 




dz 


) ] dz 


Erom above eqxiation (2.17), a generalized form of Epynolds 
equation derived by Bowson [1962} can be deduced by considering 
= 0, and = h. Bae various forms of equations of fluid-film 
lubrication derived by Reymlds [1886] , Gope [1949] , and 
ZieiSciewicz [1957] can also be deduced from tiie generalized equation 
(2.15). 


In the following seme particular cases are discuss©!, 

( l) GAS HJERICATION : KBEOSIAIHC BEAEUTG 

(ll) MJITIPIB VISCOUS lUCOMBHESSIBIE lAIHlS nJBRICAHOlT : 
SIIBEE BEAEIITG 
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(ill) THREE LiTERS nJBRIGATIOR ; HSmROSTAEEC BEASIRff 

(IV) SIEGIE YISCOJS INOOMEffiSSIBEE LiSTER IIJBRICATIOir WITH POROSITT. 

2.3 GAS nJBRIGATIOW • HIDROSTATIC BEARING 

In the case of a single layer compressible lubrication and 
considering Uie conditions 

i£.==ie.=ilL-n -5^ \ / N 

8z 9 z ~ ° P=p(x,y), n= n(x,y) 


(2. 19) 


N = ^ ’ ^2 = . 


the genaralieed equation with slip ( 2 . 15 ) i-aduoes to the following 
form 


a 3P 51 aP 


“ "*2 (">23 ' H, t I; (pu)^ + i- (pw)^J (2.20) 


where 


^o = “1 + «2 




_ (H^ H?) 

J, =»1 H, ^°2^ 




2 3 n P_ '' 2 rj 


I'w 

I 

'o 


^3 = 


( 2 . 21 ) 


(PU), =P«,(H, ptU, 1_ 
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(Pu)^ = * '■£°2 - “2 -^ -~ 

o o 

Cpt), -P«,(H, -_)_+p[v, % i 


(pv)2 - -PC^(H2 - J,^) gy +P[V2 “ “2 P 


using equations (2.21 ) and remembering "tiiat in the case of two- 
dimensional slider bearing, 


a 3H 3H 3H 3H ■ 

tPW^ = (pu)^ — ^ (py)^ — - (pu), ^ ^ . py 

(2.22) 

where Y-squeeze velocity in the (-ve) z direction, the equation (2.20) 
can be simplified as follows : 


3x n '' 3 


^ - PP,^ ap 

o 


■•■l^ H?) --^ !r ] 


2 2 r r 3y 


3z^ ^^2 “ip ^ ” PU|(H^ - J,^ )] 


p p 

4. L_ [p V (h - ~) - py (h — i)] 

3 y • 2^ 2 P 1 '‘1 P 

0 0 


- PY 


(2.23) 


Equation (2.23) can be further simplified again by considering 


= 0 , = h 
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X^ = A2=X i.e. — 


= U , = 7^ = 0 


as follows : 


(2.24) 


J.r£^ (1 + 6 1) !£i r Ei£ . 6 Xx 

as: I9n '' ° 3x ^ 37 ^ 12 n ^ h ' By ^ 


ax ‘ I2ri 


= f|^(Ph)-PV 


(2.25) 


This equation is same as Ihat obtained by Burgdorfer (1959) vdien 7 = 0 
Assuming steady state and isothenmal condition in the gas film 


we can write. 


a Pa 


( 2 . 26 ) 


Ihrther, frcm equation (2.2) we can also have. 


n 

a ’"a P 


(. 2 . 27 ) 


where suffix *a' denotes the respective reference quantities. 

UeiJog equations (2.26) and (2.27), the equation (2.25) can. 
finally be written as follows : 




X„ P„ _ ^ 3P 

1 + 6 H) _ j, 

n„ ^ P^ sx 


(1 +6-^11) 

12n h p'' 37 o 


3y 


= 1“ (ph) •- pv 


(2.28) 



mhis is Uie generalized form of Reynolds equation mth slip and 
viscosi'ty variation. 

Ifow, let us consider the flow of a gas, with constant viscosity, 
in the case of hydrostatic bearing as shown in the Pig. (2.2). Eie 
R^nolds equation for this case can be written, from equation (2.28), as 
follows ; 


5 

3 , rfa^ 

ir‘-i2n 


, ^ ® >a Pa, 3P , 

— h— )-57l 


= 0 


(2.29) 


where Pg_» may be taken as the reference values in the recess of the 
bearing 0 ^r 


iChe flow flux can be written as 
irr h' 


Q = - 


^ 6 ^a. Pa SP 


6 n 


[1 ^ 
h p 


3r 


1 


( 2 . 30 ) 


.which is a function of r. Ckxabining equations (2.29) and (2.30), it 
can be noted that 


9r 2ir ^ ^ * (2.31 ) 

which implies that in the case of compressible lubiication, 

P Q = p^ = constant (2.32) 

where Q is the inlet flux in the recess 0 < r < r. . 

How from equations (2.30) and (2-32), we have 






sigiifig;® 


pi^i 


Ml 

»c:1?^- 




— IM 
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Bie bouMaiy condition for p are 

P = Pa r = r^ 
p = 0 at r = r 

Integrating eqixation (2.33) and using the abOTm boundaay condition 
the expression for inlet flux and the film pressure can be obtained 
as follows : 

12X 


■trh^p^ 




(1 




(2.34) 


12ri ila 


and 


P 6 X ^ Jin — ^ 6 X ^ r . "T 6 X 

7 - ' — f- - ^ ^ (2.35) 

^ a * e 




^i 


Jln 


^i 


It can be seen from equation (2.35) that as X^ increases the inlet 
fl^lx Q increases for a given inlet pressure p , 


How from equataon (2.3§), we have 

r 


6X 




(!-) = : « [1 - 
Pa >“ 


in 


^i 


X„ In^ 6X ■’ 

{(1 + 12 ^)' 




In 


Jin 


e - 


(2.36) 


Since, 


< 1,.it can be noted that the expression inside the 


Jln 


square bracket of the equation (2.36 ) , is positive. Hence tT” ( — ") 

■ d p 

' ' ■ ^ . . , a :^a 
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is 3hus, it c&a. be eoncliHiCd that the pressure decreases 

as X^, the mean tree palh increases, Rirther, as 1h.e load capacity 
is the integral of the pressure, it also decreases as X^^ increases. 

2,4 miTIllLE YI scons ICTOCagRESSIBIE MYERS SJBHICAnOU t SLIDER BEARING 

It has already been pointed out that the viscosi"^ varies across 
the film-thicjkness due to thermal changes and such a study has been 
conducted by Qyale and Wiltshire [1972] . Ihe viscosity of the 
lubricant also changes with concentration of the additive vfaich may 
vary across the film- thickness, Shukla [1972] and Isa [1974] , lUrther 
the viscosity of the lubricant near the solid surface msy Increase 
considerably the closer the boundary is approached as pointed out by 
Pevenport [1973, p 19 J . 

Keeping this in view, a new concept of multiple layers 
lubrication is proposed here so that the effects of viscosity variation 
across the film and near the solid surface can be taken into account, 
Qvale and Wiltshire [1972 ] , Askwith et. al. [1966] . 3he R^nolds 
equation presented here can also be used to study the effects of 
surfactant at the bearing surface (such as fatty acids) infaich mi^t 
change the viscosity near the surface, Askwith et. al. [19^6] and 
Fein and Kreuz [1966] . lb simplify the matterjthe study here is 
confined to only multiple inccmpressible l^ers, however there is no 
difficulty in deriving a corresponding Reynolds equation for multiple 
compressible layers. 
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In tile case of ineompiessible lubricant, all G-fUnctions vanish 
^id geneialized equation with slip (2.15) siii: 5 )lies to the following forni 

a 3P g 3P 

4-tP — ] + ~ r — 1 

3x '’ 2 3x'‘ 3y ‘■ 2 ay* 


where 


®2 f 3x ^^^2 3y ^^^2 ^ ^ Bx ^’^^l 3y ^^^1 ^ 

-3 '3 (^ 2 -^ 1 ^^: 


H 


3x 


P 


3 -h i -- 4^, ' ' ^ 1 + [w]'/ (2.37) 


3y 


pf 

o 


H. 


*■0 = “i -"“2 *i 

11. 


2 dz 

n 


^0 = ^1"^ ®2 ^ 


H. 


H, 


2 dz 

n 


1 


.H. 


r, = ^ H, H. Cj H, + n 


Ij = 61 + 6^ + / 


H. 


2 z dz 


H. 


(2.38) 


H, 


P. 


p _f2z7 1\ 

2 i TrCz-T-) 


dz 


H. 


1 


^0 . ^ 1 ' 

^ = ! 7(^-iT)d2 


H. 


P, » / 


2 z dz 


H. 


1 


To study a particular case of liie equation (2.38), let us 
consider the flow of three fluid l^ers in the case of a slider 
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bearing -vrfaose configuration is shown in Pig. (2.5). Considering also 
the effects of slip at the surfaces in this case asd roaeaibering eq.xiations 
(2.1) and (2.3)j we have 

6 

where is the viscosity of the lubricant at the surface and 0 is the 
coefficient of sliding friction at liie surface. It is noted here, that 
as 0 increases X decreases and 0 tends to infinity in the case of no- 
slip at Uie surface. How, keeping in view the' plqrsieal situation 
represented in Pig. (2.3) and considering the following 


= U, Ug = = Tg = 0 


“1 ' “2 = ®1 = ®2 = i 


(2.39) 


P = const j 8^ = hg = const 

the E^nolds equation applicahle to this situation can be deilved 
from equation (2.37), as follows : 

a ap a 
S- fp — 1 + r p - — ] 
ax ‘ 2 3 x ^ ay ^ 2 ay 


where 


- (u)j ] 

+ D |- (r^) + [m]**'® 

ax Pq •• ■'o 


O 112 ^1 ^ 


(2,40) 
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■B ^ ^ ^ 

2 12^2 12n^ 


^ /h+Hv /h. H V 

^3 = ^— > 


(u) = _i (MI)!£ + _s_ 

vu^g B 2 ^ as BP 


^2 = -r(¥)F 


(2.41) 


= (u)^ g- (h+H) + (t)^ ~ (h4fl) - 7 
SiEq)lifyirig further equation ( 2 ^ 40 ) by using equation (2.41 ), we get. 


3 3P 3 3P U 3 

3x f ^4 "ix ^ ^ 3y ^ ^4 3y^ ~ ¥ 3x 


(2.42) 


where 


P. 


4 i2n 


h^ H^+3H^+3Hh^ (IHH)" 
f2Ti„ 12 TL 2B 


''2 ' 2 Tij 2 B 

Applying equation (2.42) to the case of one dimensional beaidng, we 
have ■ 


(2.43) 


±_ 

dx 


'*4^’ 'II; (h«)- V- 


(2.44) 


Integrating above equation (2.44) with the boundaiy condition (keeping H 
constant) 


we get, 


dp 

rr- = 0 at X = x„ , h = h 
dx o o 


dp -^(h-h^) + v!V(Xq-x) 
dx P. 


(2.45) 


(2.46) 
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following the usual procedure j -foe load capacity and the 
force of friction in this case can be written as follows s 


dp 


U. 


W = / (-X dx = - / 


0 


(h-h^) + 

~Tr " 


dz 


(2.47) 


^ = / %(B 


^3z'2=0 


dz 


= / 


Ns 

0 ^o 


h-ffl\ 2. 


(^) 


+ V(x^-x) 


“^4 


3 dz 


(2.48) 


Prom the expressions of P^ and P^ it can be seen that P^ and P^ decrease 
as 3 f or T]^ increase. Hence from equations (2.47) and (2.48) it can 
be concluded that the load carrying capacity and the frictional force 
increase as 8, n^or ng increase. 

In the case of pure squeezing between two parallel plate, equation 
(2.44) can be integrated with ustial boundary conditions to give the 
expression for film pressure and load capacity as follows ; 


P =2^ x( lr-x) 
L 

W = / p dx = 
0 



(2.49) 


(2.50) 


iis in the previous case it can also be seen that load • ; 

carrying capacity -increases a§ 8, Hence, it is 

concluded that the load capacity and friction force increase as the 
viscosities of the layers as well as the coefficient of the sliding 
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friction increase, lElie effect of slip is to decrease the load capacily 
and the force of friction, 

2,5 BIREE lilERS PJBBICAIIOK t HYBROSIAiEEC BEJEMG 

Consider the flow of three incompressible fluid layers in the 
case of externally pressurized hydrostatic bearing as shown in Eig. (2.4). 
Bie equation governing the flow in this case can be calculated as 
follows ; 


^ g] =0 (2.51) 

where is defined in the equation (2.45). 

Integrating equation (2,51 ) with respect to r and using the 
boundary conditions 


we have^ 


P = P. at r = r. 
1 i 


p = 0 at r = r 


0 2? 

p _ p .. £,n 

^ 2irE, r. 

4 2. 

p. = ' • 

2irP^ r^ 


aSie load capacity can be evaluated from the expression 


(2.52) 


(2.53) 


2 ^e 

W = iip^ r^ + J 2 ir r pdr 


(2.54) 


as follows : 


Q(r^ - 4) 

4ir P, 


W = 


( 2 . 55 - 2 ^ 











or 


/ 2 2x 

~ 

2 iln(i) 


(2.55-b) 


It has already been noted from the equation (2.43) that decreases 
as 0 , the ^efficient sliding friction increases. Hence, it can be 
seen from liie equation (2.55-a) that for a fixed flux Q in a hydrostatic 
bearing, the load capacity increases as viscosities of the layers increase 
and it decreases with slip at the surface, 

2 £ SIHQ-I.E 7ISCQUS IHGOMPRESSIBLE LAYER lUBHiaiTIOH WIB I POROSIff 

Let us now consider the case of single incompressible layer 
lubrication with slip velocily at the porous boundaiy. Taking, 


= 0 » ^ ~ '^1 ~ ^2 ~ ~ const- , p = const. 

Ij = 6.J = 0, Xg = ^2 = X 

X 

i.e. a/= 8- - 0, a_ = 8^ = *- 


(2.56) 


'1 "1 
H, 


n 


the EeynoldS' equation (2.57), with slip in the two dimensional form 
can be written as follows : 


I p — ] + .i_ r p — "i 

^ 2 3x 3y ^ 2 3y r 


3x 


Ch(u)^ 


9y 


[ b(v) 1 


(U - U. ) E 

3 r 2__! 3] ,7 


f-t 

ax 


s 


(2.57) 


where 



45 


(u), = (h -— ) 5- + [tJ, -i (-21)] 


F' ax 

o 

o 

J, J, ^ 

o n ’1 2 ti 


5- = id- ' _ ^ 

2 I2n + X ^ * 3 2n 

<1, ap* 

''s = %- <35-U ’ 


and p* is the pressure in the porous region. 


Sin^ilifying equation (2.57) by using expressions (2.58) 


have 


3 ap 


9 r h 


3 ap 


[f--] E-rr 1 

2 


ax ax ay ay 


^ a r \ TT , h(h + 2X) „ 1 

i--(hTTr 


- 12 » 


where 


•i-- tl +-Ar ] 


1 


1 + 


In this case, on comparing equations (2,l) and (2.4), 
1; ■ 
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X = f and the equation (2,60) can be rewritten as 


X 

0 , 


= {1 + 


[1 + 1 


} 


( 4 )' 


(2.58) 


we 


(2. 59) 


(2.8O) 


(2.60-a) 
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It is noted here that when the slip coefficioat S "*■ 0 whxle 

when 5 •*■ “ » "*■ 1 (ifo Slip), 

Now, it can be noted 1±iat the equation (2,59) reduces to the 
case of two dimensional squeeze film between porous and a solid 
rectangular plates as discussed by Wi [1972] . 

In idle case of eicteinally pressurized porous squeeze film 
bearing, Wang [1975] , 

n '‘82 ^z=h n ^ H 


where H is the thickness of porqus facing and p is ihe externally 
applied pressure, the Reynolds equation (2.59) can be modified to. 


2 2 
8 P 8 P 

— 5- +— 2 
8X: 8y 




n 



(2.61) 


In the following, we study the case of oi^-dimensional erternally 
pressurized porous squeeze film as shown in Pig. (2,5). 



» 
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Solving equation (2.65) by using boundary conditions 


p = 0 at X = 0 


I. 

■S * 0 -at X = - 


(2.64) 


we have. 


P = [P^ " 


cosh e(z - —) 


S , 2 ' 


■ e ^2 • , KI 

K cosh ~ 


1 


( 2 . 65 ) 


fEhe load capacily of liie bearing is foxind by integrating the piessanre 
in the film as follows r 

K^L 

(2.66) 


L g tahh 

W = / p dx = I [pg [1 - - ] 

K! 


In dimensiohless form it yields 


EL 

2 


W = 

where 


^ <}> 


( 2 . 67 ) 


£ =1- V w= ^ 


^Pe ' 


2 -^ 


. dt ni 

5 ® , 


v = 


H h"^ 


(2.68) 


GSae time of squeezing for a constant load can be found as follows : 


5=1 / 


dh 


( 2 . 69 ) 


¥ 


h [ 
1 


- 1 ] 


where h^ “ h" 


V 
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When o ^ =1 ( 5 =«o -tije expressions for ¥ and I reduce to the case 
corresponding to no-slip condition. !Btie load capacity and the time of 
squeezing given by equations (2.6?) and (2.69) are plotted in KLgs. ( 2 . 6 ) 
to ( 2 . 11 ) for various parameters. Heie liie load capacity and the 
time of squeezing are calculated for 5 = 0 . 1 , 5 = 0.5 and 5 = ® . 


It is clear from ELg. ( 2 . 6 ) that for no squeezing (y = o) the 

load capacily W decreases as — inc^reases aM increases as §(the slip 

coefficient ) increases for a fixed 4 • Sit in the case of squeezing 

(F *s 3»0), Esee ELg.(2.7) 3 } W decreases for a certain SEisO.1 values of 
H H 

— and -Hiereafter increases as ■“ increases. lUrther it is also noted 
L Ii 

H 

from Fig. ( 2 . 7 ) that for small values of the load capaoily is lower 
for small values of I” vdiile "this situation is reversed for hi^er values 
of U • Frcmi Fig. ( 2 . 8 ) it can be seen that for -j- = 0.05, the load 
capacity increases as 41 increases for Y = 0 but it decreases for V = 3*0- 


It is observed from Fig. (2.9), ( 2 . lO) and ( 2 .II) that the 

existance of slip at Ihe surface shortens the time of squeezing between 

H 

two plates. iOr a fixed y ~ if con be seen from lig, (2.9) that 

N - 

the time of squeezing increases as ■g— and decrease . Further, 

o 

from Fig. (2.9), it is noted that the time of squeezing increases as 5 
increases. In Fig. (2. 10), it is shown that the time of squeezing 


H 

I 


increases as -r increases for a fixed -r^ = 0.25. From Fig. (2.11 ) for 


1 


a fixed = 0.25, it can be seen that the time of squeezing increases 

■ o ■- ■ ' . H ' ' ' ' ' 

as <|) decreases but the value corresponds to 'r ~ i® still hi^er 

Jj 

than that of ^ = 0.01. 
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2.? RESULTS 

In this chapter, a generalized form of Reynolds eq^uation 
applicable to fluid (gas as well as liqiiid) film, lubrication is 
derived by considering the variation of fluid properties across as 
well as along the film thickness with slip velocities at the bearing 
surfaces. Various particular forms are discussed. 

In case (l), a Reynolds equations with slip is derived fcxr 
compressible gas lubricant and the case of an externally pressurised 
bearing is discussed analytically. It has been shown that the pressure 
in the film decreases as the mean free path incareases. 

In Case (ll), three different layers of fluid with different 
viscosities have been considered and a Reynolds equation wilh. slip for 
two dimensional slider bearing is deduced. In one dimensional case, 
it has been shown that the load capacity and frictional force increase 
as the slip coefficient and the viscosities of the fluid layers 
increase. The same result has been pointed out in the case of an 
externally pressurised bearing. C£x.,ie. ( li 

In general it has been noted that the effect of slip at the 
bearing surfaces is to decrease the load capacity and frictional force. 

In case (l^), the slip velocities are considered at the porous 
surfaces to obtain ageneralized Reynolds equation for incompressible 
lubricant. Even, in the case of externally pressurized porous squeeze 
bearing, the existant of slip velocity at the porous surface reduces 
■fee load oapacity and shorten the time of squeezing. 
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3c,y,z 

p 

X ,X , S ,5 - 
1 ’ 2 * 1 ’ 2 

( )2 


LIS I OF IiI03!Aa?I0HS 

film thicloiess 
initial film thickness 
final film-thickness 

distance between liae surfaces and the origion. 

length of bearing 

p3?essure 

inlet pressure 

external pressure 

radial coordinate 

inlet radious 

outlet radious 

time of squeezing 

velocity components of the film in x, y and z directions 

velocities of the surfaces at z = Hi and z = , in 

X direction 

velocities of the surface at z = Hi and z = 
y direction 

Squeeze velocity 

load capacity 

rectar^ilar coordinate 

density of the lubiicant 

viscosities of the lubricant 

parameters defined after equation (2.8) 

pomeable parameter of porous facing 

value at z = Hi 

value of z = Hg 
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CHAPTER - III 


A THSOSI OE STOCHASTIC IHBRICAHON 

5.1 STTROMOTIPIT ; 

Owing to machining limitations, it is not possible to 
manufacture perfectly smooth surfaces of the bearings, naturally, 
therefore the question arises as to what happens to the vailous 
characteristics of the bearing when the surfaces are roug^. This 
question has been taken up by some workers Burton [1963] , Shukla 
and Prasad [ 1966 ] , Dowson and I'/homes [ 1971 ] in the past assuming 
that rou^ surface to be represented by a single sine (cosine) 
wave or a series of sine (cosine) waves and thus modifying the 
film thickness in the usual study of the bearing characteristics. 

It has been pointed out that the load capacity, frictional force 
etc. can be different from their value corresponding to the aaooth 
surface and this difference depends mainly on the amplitudes and 
the wave lengths of the waves representing the rou^ surface. !Siis 
procedure, called the deterministic approach has also been applied 
to study the characteristics of rollers Dowson and Whomes [ 1971] 
and spiral groove bearings Wildmann [l968] . 

In another method, called the stochastic approach, the 
surface rou^mess is assumed to be represented by a stochastic process 
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and the usual procedure of study is followed by taking the statistical 
mean of the basic equation. In fact, it has been daronstrated that 
the density distribution of rou^ness heists in usually machined 
surfaces follow nearly the gaussian distribution and the gap between 
two rou^ surfaces can be represented by a stochastic process 
Papoulis [1965 ] . Using this approach the characteristics of an 
infinite slider a.nd short journal bearings has been studied Szeng 
and Saibel [ I967a, 1967b] , Bae Reynolds equation applicable to 
finite rou^ bearings has also been derived under certain assumptions, 
Christensen [ 1969-70] , Christensen and Tender [ 1969a ] and the case of 
finite slider bearing has been investigated, Christensen and Ponder 
[ 1970] . Biis concept has also been extended to rough bearings using 
compressible lubricants, ELord [ 1973 ]. 

In this chapter, we derive a generalized foim of Reynolds 
equation applicable to finite rcxigh bearings by using stochastic 
approach, Sveshnikov [l966] and lin [1967] Bie cases of hydro- 
dynamic one dimensional step slider bearing and hydrostatic bearing 
are investigated. 

5.2 BASIC IBUATIORS 

Consider tiae flow of an incompressible lubricant between 
a stationajy rou^ surface and a moving smooth surface (see Pig, 3.l). 


Th e equation governing the pressure in the film is given by. 




(3.1) 


a 

ax 




ax^ 


4* 


ay 


= enU-B + SnTifi 


-M 

ax 


an 


where the piTessure p can, In principle j be determined for a given 
film thickness H. Since the film thickness function H is a 
stochastic process, the pressure also becomes a stochastic process 


and can only be determined by an averaging process, Pefclenik [1968 ] 
and Papoulis [ 1965 ] . !Ehus, for a given random input H, equation (3.1 ) 
is a stochastic differential equation for deteimining the random 
function P, Christensen [ 1969-70 ] . 


How, by taking the statistical average or expected value of 
each terms of equation ( 3 . 1 ) we have , 

h +37 

where 

00 

E(s) = / sP(s) ds and P(s) is the probability density 

distribution of the stochastic variable s, H being the film- 
thickness function, Tseng and Saibel [1967 -a] and Chris tens|eai[ 1969-70] 
written as follows : 

H = h(x,y) + hg(x,y,6) (3^3) 

Here h is the nominal film-thickness (the film-thickness 
corresponding to smooth surfaces) and h^ is the part of the f ilm - 
thickhess due to surface rou^ness measured frcm the ruominal lev^ 
and is assumed to be a (stationaiy ) ergodic, stochastic process. 
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In general, two stochastic processes P and H are correlated 
and even if all the statistical chaiacteristics of H are known, it is 
not possible to find out the inean of such terms as ~ in terms of 

oX 

3 P 

the means of H and — explicitly unless some extra relations are 

o X _ 

assumed, Christensen [1969-703 , Christensen and Otonder [ 1969a, I969b ]. 
As such, since a stochastic function can be approximated as a series 
expension in terms of deterministic and dependent or independent 
stochastic functions Sveshnikov [1966] and lin [1967] , it is assumed 
that the flow fluxes which are stochastic processes can be written in 
terms of H as follows : 


^x 2 ^ 12n ax 


ap u 


M 


V 


1 ap 


-f I 

i=0 




Q = I B. 

7 2 I2n ay 2 i 


(3.4-a) 

(3.4-b) 


H 


ere A^*s and ’s are non-random functions of x,y,U and V or 
constants depending upon the lype of rou^ness, rou^ness slopes, 
minimum film-thickness etc. and may be determined Iheoretically or 
experimentaly [ see appendix lY ] . However , it wjl?- -SQaa latpr- that 
by choosing A.- 's and B. ’s suitably, some well known cases can be 
obtained, Tzeng and Saibel [1967-a] , Christensen [ 1969-70 ] , 
Christensen and Ponder [l969-b,l970 ] . 


In general, as pointed out by Ponder and Christensen [ 1972 ] , 

the flux at a point in a lubricant film may not be determined by the 

local film- thidmess alone such as in the case of waviness (long wave 

undulations ) where h is comparable to the nominal film-thickness. 

s 
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However, in the case of roughness i. e* the case of hi^ frequency 

undulations with smaller amplitudes (i*e* h « h) the rep3?esentatior]s 

s 

(3*4) are applicable. Parther, it may also assume that the mean of 
hg is zero [i.e. E(hg) = 0] and then E(h) = h. 

IGaking the expected value of equation (3*4-a) we have 


U 


e(h) — n e(H^ — )■ = 2. 
^ ^ 12n ^ 2 


I e(H^) 

i=0 


(3.5) 


Again, dividing equation (3.4-a) by IT and taking the expected value, 
we get. 


U 

2 


E(^) - 

H 


ik = 


M 

I 

i=0 


A. 

X 


El 


(h^-5) 


(3.6) 


Eliminating where 0 ^m from equations (3*5 ) and (3.6), 

we have 


U 

2 


" 2 ^ 1 A^E(H^) + . 


i=0,i?4a e(^'”^) 




. n-i L 

ET i=0,i^m 

Similarly, from equation (3»4-b) we get, 


(3.7-a) 




^ E(H) - -rii- E(h5 |~) = I [ I B, e(h3 ) + 

3=0,j7^n 


2“'“^ 12n 3y' ^3 ^ ■ e(h^"^) 


if " e(b“ » 

(3.7-b) 



Uow from equations (3.2), (3.7-a) and (3.7-b) the generalized form of 
Reynolds equation, applicable to a finite- bearing of 'vihich the 
moving surface is smooth, can be written as follows : 


3X ^ 'J sy ‘ ay ■'' 


= enu 


3 z 


[ I A E(H^) {Efe) - I A E(H^“^)}] 


i=0,i74i 


+ 6nv|--[ I B E(H^) + .. X^) {e(-^) - I B E(H^“^)}1 

(3.8) 

Einally, equation (3«8) can be written in a more simplified form thus ; 


-|- [ h^ f |- e(p)] + 1-[ h^ f e(b) ] 
ax'- m 9x ^ ^ 3y n 3y ^ 


-s K'> + «™9T 


(3.8-A) 


where 


f 


E(lf) 


E(if5) 


M 


I A,V)+v3{e(-L) 

i=0,i?4ii H 


M 


i=0,i^ 


and similar expressions for f^ and g^. She functions f^, f^, and 
g^j would have different functional relation wi-Ui e(H^) for different 
types of rou^ness and indices m,n depend upon rou^ness heists, 
rou^ness slopes etc. ihe solu-tion of equa-tion (3.8-A) can be 
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easily plotted for different values of these functions in tabular 
form for convenient use. 


iUrther, in the case of compressible lubricant, if p is 
assumed to be function of the mean pressure e(p) only i.e. P= p{e(p)} 
and its variance is neglected, then by follo?.dng the above procedure, 
the corresponding Reynolds equation for gas lubricated bearing can be 
deduced in the first ^proximation as follows ; 


3 r 


P{E(T)}|^E(P)] [ 


e(h''~^) 


PtE(p)>|^ 


E(P)1 


= 6hu|-[ I JLE(H^) {sfe) P{E(P)} ~ I A E(H^*^ )}] 

"" i=0,iA E(lf If i^O,i^'" 

+ 6r)7|-[ I B E(H^) {e(-^>{e(p)1> I B.E(h^"^)}] 

3=0, ^ 3=0,3,fe J 

(3.8-B) 


Here A; 's and B. 's represent the various "type of rou^ness,rou^ness 
distribution and rou^ness correlation of the bearing surface. 

How in the following some particular cases are considered by 

taking B. = 0 for all j and V = 0. 

J 

3.3 PARfICTJIAR OASES 

CASE (i)* ^ considering A^ = 0 for all i in equation (3-8), we 

get the E^nolds 1ype equation applicable to correlated r outness 
along sliding direction as follows : 
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S r E(#) 



1 3 r E(H^) 3 

37 E(H-^-3)37 


e(p) } 




(i-l) If m =0 and n = 3 , we have from e'-uation (3.8-c) 


(3.8-0) 


e(^) 


•h=' 


( 3 , 8 -D) 


which is a R^nolds type equation for non-uniform, transverse rou^ness, 
correlated along sliding direction (one-dimensional transverse rou^ness) 
Christensen and fonder [l970 ] . 


In one-dimensional case the above equation reduces to the 
following form vshich is applicable to an infinite slider bearing, 
Tzeng and Saibel [1967 -a] . 


t-J 

dx TT^/I ^ <t3C 


E(h3) 


e(p)] = 6hU[ 


E(\) 

H 

IT 


( 3 . 8 -E) 


(i- 2 ) If m = 3 and n = 0 , then ffom equation (3.8-c), 


3 r 1 3 


A[E{H5)nE(p)] 


e(p)] 


= 6nU |- [e(h^) e(J^) 1 

3x g2 


(5.8-E) 



which is applicable for non-uniform, longitiidinal rou^ness, correlated 
along sliding direction. 


(i-3) If m = 3 and n = 3» then from equation ( 3 . 8 -c), we have 

[e(h^) s(p)] +|~ [e(h^) ~ e(p )3 = 6 nu [e(h^) e(^) ] 

(3.8-ff) 

which is a Reynolds type equation for uniform, longitudinal rou^ness, 
correlated along sliding direction. 

(i- 4 ) If m= 0 and n = 0 , then from eouation (3.8-c), 


8 X 


E(^) 


3(P)]+|r[ 


1 


3y E(l^)3y 


e(~) 


] (3.8-H) 


which is a R^nolds type equation for uniform, transverse rou^ness, 
correlated along sliding direction. 


CASE (ii) ; considering =1, A^ = 0 for all i in eqxiation (3.8), 

Vi ,m 

we get the following Reynolds type equation applicable to less.- 
corarelated rou^iness along sliding direction. 


8x ^ ^ ay j5(H^“^)3y 




■ 3 \ 3x 


ax 


(3.8-1) 


(ii-l) : If m = 0 and n = 3» ftom above equation, one gets 


:[ 


1 a 


ax'- \ ax 


E(p)] + 1 ^ [ E(#)|^ E(p )3 = [e(h) ] (3.8-j) 


ay 


ay 



71 


which is applicable fcxc noct-uniform, ■fcran.svexse rou^iness, uncorrelaied 
along sliding direction. 

(ii-2) : If m = 3 and n = 0, from, equation (3.8-1 ), we get 

^[Kh3)|jE(p)] +|_[-^iKP)] = 6w|j[E(H)] (3.wt) 

which is applicable for, non-uniform, longitudinal rou^ness, un- 
correlated along sliding direotion (one-dimensional longitudinal 
rou^ness), Christensen and fonder [1970 ] . 

(ii-3) : If m = 3 and n = 3, from equation (3 *8-1), one gets 

~ [ E(h^ ) 1^ E(P) 1 + Ij: [ E(H^ ) ^ E(P)3 = 6nU 1^ [ E(H)-J(3.8-I) 

which is applicable to uniform, longitudinal rou^iness, un-correlated 
along sliding direction (isotropic) as derived by Christensen and 
fonder [ 1969-a] . 

(ii-4) : If m = 0 and n = 0, from equation (3.8-1 ), we get 





Mp)] + 



which is applicable for uniform, 
along sliding direction. 


—^—2(^)3 = 6nu |-C e(h)3 (3-8-m) 

) ay ax 

transverse rou^ness, un-correlated 


It is seen from the above forms of equations that when m and 
n are equal, the term 'uniform rou^ness' has been used, ^ile when, 
they are different ,the term 'non-uniform* is suggested. The term 



‘transverse roughness’ basically implies a correlation between the 
film- thickness function and Hie pressure gradient in the x— diiection 
while the teim longitudinal rougjiness suggests non-correlation in the 
above terms. Ihe correlation or otherwise between different film- 
thickness functions in the sliding velocity term has been indicated 
in various cases by writing these relations in the "sliding direction". 

It can also be noted that tte rou^ness correlation decreases 
as m and n deviate from zero to three and as increases from zero 
to unity. 

In the following, the two particular cases of bearings are 
studied where norminal film-thickness is constaat. 

(I) S 2 EP BEARDJ& 

(II) HIDROSIATIC BEABIN& 

3.4 SfEB BERING 

Consider an infinite rou^ step bearing as shown in Eig. (5.2). 
Since, for all practical purpose e(e^) and E(i^) can be assumed 
constants, we have from equation (3.8) t:he equation determining the 
pressure as follows : 

A 

5 — E(p. ) = i = 1 in region I 

o X X 


i = 2 in region II 

Since the pressure at the step is continuous, we have. 


(3.9) 
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Using the continuitjr of flew fluxes in the two regions and 
remembering equation (3. 7 -a) we can write : 




E(if) 


E(l^~0 


M 


i=05i74i 
i-3, 


E(i^"^) 


{E(^)- I A.E(Ht"^)}] 


E(I?) 


I2n 




Ub 


M 




i= 0 , l/m 


E(l^) 


M 


{eS) - i A^E(Hj“^)} ] 

i= 0 ,i/ 4 ii 


(3.11) 


SolTTing equations (3.IO) and (3.II) and integrating the resulting 
expressions under usual boundary conditions, the following equations 
for determining the pressure is obtained. 


and 


E(P^ ) = 


eCp^) = 


6Un (i-k) (g^ - gg) X 
(l-k) f^h^ + kf^ h^ 

6Un k (g^ - §2^ (E-x) 
(l-k) f^h^ + ^ ^2 4 


0 < X < kL 


(3.12) 


kL < X < L 


(3.13) 


where 






M . , , M 

= . J ^ - . I. , A .2 

i=0,i;4i Hjj i=0,i?4a 


(3.14) 



Proceeding in the usual manner, we write the following 


expressions for the flow flux and the load capacity respectively 


as follow 


E(Q) 


(i-iOCg^-g^) 


'5.15) 


(l-lr)+ 


E(w) = 3hUbl' 


2 k(l-k) (g^ - g^) 


(l-k)f^h^ +I':f2h2 


The frictional force on the moving surface can be written 
similarly as 

3k(l-Tk)(g -g )(f’h - f 'h ) - 

e(p) =nbUI[ +{kg’ + (l-k)g’}] (5.17) 

(l-k)f h^ +kf ^ . 2 


where 




. = E(^). 3 [E(n). f 


=a 


i=0,i74a 


- - I. , ^1 . “= 1.2 


i=0,i;4i 


(3.18) 


The expressions for f^ , fj , and g^ are expressed 'in 
terms of a/h^ for various values of m and A. ‘s in Appendix I. 


The friction-coefficient can ^c defined as follows : 


"If} 


(3.1°' 



The expression for e(w), eCe) and e(v) are plotted in 
T^gs. ( 3 . 3 ), (3.4) and ( 3 . 5 ) respectirely for the following 

N 

combinations of m and AJs (-r; — = 1*656, 1: = 0.7 ) for different values 

2 

of cr/hg- 

(a) JLj = 0 

or for n = 0, A. =0 for all i 

i 

= 1 i/'O 

(b) A,=o-l 

' i 

or for m = 3j = 0 for all i. 

Aj = 1 i/^3 

Hence, m = 0 corresponds to the transverse rougjiness while 
m = 3 corresponds to the longitudinal rou^ness. The condition 
A*| = 0 implies a correlation between the film-thidcness functions 
in the sliding velocily tem while A^ = 1 corresponds a non- 
correlation in these functions. Physical3.y, for the same asperity 
hei^t, the case A.| = 0 corresponds to sharp and rough asperity 
slope with small asperity base T/hile A* = 1 corresponds to wavy, 
less sharp, smooth asperiiy slope with comparatively large base. 

Erom the figures (3.5)> (3.4) and (3.5) it can be seen that 
in case (a) the load capacity, friction force increase but the 
coefficient of friction decreases as o/h^ increases. It can also 
be seen that the increase in load capaeily is mo 3 ?e significant than 
that of friction force and this leads to a lovrer coefficient of 


friction, ffiiese results axe similar to that obtained by 









Tzeng and Saibel [ 1967~a3 for rougji slider bearing (i.e m = 

=0, Ajj^ = 0 for all i). It can also be seen that the load capacil^" 

1^0 

in friction force are greater in the case of a step slider bearing 
than the inclined slider case, for the same rou^ness gmd other 
characteristics, Tseng and Saibel [ 1967-a] . 

itirther vhen A^ = 0, the increase in load capacity of the step 
slider bearing in comparison to smooth step bearing is about 35^ for 
o/hg = 0.3 while the decrease in coefficient of friction is about 4^. 

To see the effect of A^ in case (a) it is noted from Pigs. (3 * 3)5 
(3.4) (3*'5) that the load capaciiy, friction force and the 

coefficient of friction are greater in the case Aj = 0 as ccmpared to 
the case A^ = 1 for given 

In case (b), these results are opposite to that of case (a) 
in all cases except the case A^ = 1 ,m = 3 > where ihe change in frir+'’ o • 
force is not appreciable ?ath respect to the increase in <^/h 2 * 


3.5 HYDROSTATIC BEABIFG- 

How consider "^he case of a rough hydrostatic bearing with 
constant naninal film thickness as shown in Pig. ( 3 * 6 ). Assuming 
that the variance of the film-thickness is constant we have from 
equation ( 3 . S') the equation determining the pressure in this case as 
follows : 



E(lf ) 




(3.20) 
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Integrating equation (3.20) and using the definition of 
flow flux, 


E(Q> = 


TTx" 


E(lf ) 3 


e(p) 


(3.21) 


E(if 5) 

Solving ecjuation (3.21 ) and using boundaiy co-nditions e(p) = e(p. ) 
at r = and E(p) = 0 at r = r^, the expressions for pressures 
can be written as follows : 


E(p) = 


6nE(Q) 

— =- — In — 
Tfh^f ^ 


e 


6nE(Q) 

E(p ) in^ 

irh^ f ^ 


(3.22) 


(3.23) 


where 


f = 


E(lf) 


^ ^ and X = — 

° i? ''e 


Ihe load capacity in this case is given by 

3nE(Q)rJ ■ 

- E(B) = (1 .X=)^ 

ll ^O 


( 3 . 24 ) 


!Ehe behaviour of the function is shown in the Pig. (3*7). 

■"0 

Erom this, we see that this function increases for m = 0 and decreases 
form = 3 as n increases. Phus, ws in:? er from equation (3.21 ) that 
the flow flux, at a given inlet mean pressure l(Pj^), increases for 
m = 3 and decreases for m = 0 with the increase of a . But if e(q) 
is assumed to be constant, then the mean pressure E(p) decreases 
for m = 3 and increases for m= 0 as o increases. Similarly, for 
a given mean flow flux, the load capacity decreases for m = 3 snd 
increases for m- 0. 




We have also studied the characteristic of this bearing by 

representing the rougbness as a series of cosine function 

[appendix II ] for All the characteristics obtained in this 

section for m = 0,. are qualitatively same as that of obtained in 

appendix II. But if the variance a is chosen such that 
M 

a — j then all the characteristics are exactly the same 

quantitatively as well [ Por example compare eouatioas (3.24) and 
(3.46) for ^ « 1 and the appendix I ]. 

Ihe frictional torque, #iieh is defined by 
4 

irriJJr . 

E(t) = — (1 - r) E(I) (3.25) 

-j 

also increases as 0 increases as liie function e('^) always increases 

H 

with o . 

a 1 

Einally for ■^=•2, m = 0, the fmction ~ increases by about 

^o 

30^ [see Eig. (3.7)] > hence -tie load capacity also -[see equation 
(3 *24) ] increases by the same amour t. 

3.6 DISCUSSION 

In this chapter, the generalised Eeynolds equation applicable 
to rou^ surfaces has been derived by assuming that the film-thickness 
function is a stochastic process and satisfies the basic Eeynolds 
equation. To find the average of each term of the stochastic form 
of Reynolds equation, it has been assumed that the fluxes which are 
stochastic functions can be represented by power series of the 



stoch^tic rariable H, mth variable coefficients* If*om the resulting 
equation it is seen that, since H is noimal, the mean pressure can be 
expressed only in terms of the variance o (roughness parameter) as 
the hi^er moments of H can be expressed in teimis of this variance* 

"By choosing different values of A., B., m and n ?fhich in general 

j 

depend upon the type of rcu^ness (i. e. mean, variance, asperity slope, 
asperity radious, elastic and plastic parameters of the surface etc.), 
we can derive the various particular cases of the R^nolds equation 
applicable to the 3?elevant rough bearing surface. Erom the generalized 
Reynolds equation (3*3) all the known cases have been derived as 
particular cases by choosing different values of parameter involved. Surthe? 
it has been pointed out that one ean get. the' various- forms, of ^ . 

Reynolds equations applicable to correlated and un-correlated rou^nes? 
for different values of m,n and A^'s [ see CASE (i) and CASE (ii)] . 

In the case of rough step bearing it is shown that the load 

capacity and friction force increase as o/h^ increases but the 

coefficient of friction decreases for the case m = 0, A^^ =0 

i/Q 

!Qae case of hydrostatic bearings also investigated both by 
stochastic and sinusoidal approaches [ see appendix II }. In both 
the cases it is shown that the load capacity aid the frictional torque 
increases as the rou^iness parameter increases for a constant mean 
inlet pressure, but the flow flux decreases. It is possible by a 
suitable choice of variance and wave laigtii relation, to make the 
various characteristics of this bearing exactly the same in the 
two approaches. 



APPEEDIX - I 


” APPROAli!ai‘rE DEPEElI»m-?ATIOH OP PUKGTIOIIS 


In the case of hydrodynamic lubrication, the ratio 
be assumed to be small (|- « i ) and we can write 


h 


may 


if" = [h + h j 


m 


h 

^h 


[1 (i) +3!!fcl) 


2 +... ] ( 3 . 26 ) 

Assuming the density distribution function of h to be normal with 


variance, a , one can write from equation (3,26) «1 


E(lf) ,11“ [1 4^/(0)] 


(3.27) 


h 


Phis function increases as a increases for m = -1. Using equation 

( 3 . 27 ) the function f can be approximated as follows ; 


f = 






(3.28) 


-"a 


This function increases for m = 3 and decreases for m = 0 as o 

^2 . . 
increases. The function ■tp" can also be approximated to the following 

1 

equation 


(i^? - ^-o) 


^ - [1 + 3(m-2) " I ' 2 ~ ^ ' + *^^(°) ] 


( 3 . 29 ) 


1 


^2 


Since h^ > h^ (in the case of step bearing), this function also 
increases for m = 3 s^nd decreases for m = 0 as o increases. 

Also, the. function 

— ^[-1 + (^- 3 l 5 - + a'^(o) ] (3-30) 


f ’ 
a 


h 


is equal to 1 for m = 3 and it decreases as o increases for m = 0, 



The function 


(3.31) 


f,'b, 

12 

is independent of 0 fop in. = pjivi 

xa j am It increases as 0 increases for n = 

Also as the function g is, 

M G 




we have, 


(3.32) 


^ct 


(i) m—O, A.| = 1 » ^ ~ ^ i» 

i/^ 0,1 

for 

(ii) m = 3, A, = 1 , = 0 for all i 

i;^1,3 


and. 


E(-V) 


ii^ “ 

^ ^ 2 
(g^ ~ §2'^ ^ ■ ^2^ + '^'^(0)3 

for m = 0, Aj_ = 0 for all i. 5H1US the function (g.^ 
i/0 

as a increases • 


- gg) i 


(3.33) 


increases 


Again, 


and 


= E(i^) . E(^) = tl + 6 ^ + 0'‘(o)] 


(g,-g2)» (vh) [1-1^ -“"(0)1 


(3.34) 
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for m = 5j - 0 for all 03ius the fmction (g -g 
increases. 


Siniilarly from the function g' , 


M 


i ivi • 

-E(^)+3[E(^)- I 

ot a 1=0, i;4i 


M 


- - . 1. ^(1^-5) 

IT i=0,i74i 


we have for m = 0, A^ = 0 for all i, 

g'=4E(^)-^— {A)}^ 

“ ^ 

a 


which increases as o increases. 


Prom equation (3.35), for m = 3» -^ = 0 for all 


g 


r = 4 e(^) _ 3 e(H ) E(-^) 


a 


^ / "V O'' 

“ t 

^[1 -5^.aV)l 

o h 

a 

which decreases as a increases. 

P3X>m equation (3.35), 

(i) for m = 0, A. =1, A. = 0 for all i, 

i;^,1 

and (ii) for r. = 3, A| = 1, A^ = 0 for all i, 

i^1,3 
.2 


'I.A 1 


gd = ^r'> = - 

a a 

which increases as c increases# 


[1 +2; + /(o)i 

4 


) decreases as 


}] (3.55) 


(3-36) 


i^ we get, 


(3.37) 


(3.38) 
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■^PMBIX-II s SriTUSOIDiL APPROACH FOE HYEEOSTITIC 


In the case of a hydrostatic bearing [ see Fig. (3.5)1 tife. 
equation determining the pressure is given by 


dp _ 6nQ 

dr " " 3 

■nva 


(3*39) 


which, on integrating and using the boundary condition p = 0 at 
r = r^ and P = at r = r^, give 


p = m p _n 


(3.40) 


tS 


(3.41) 


Phe equation (3*40) determines the pressure while equation (3.41 ) 
gives the flow flux in the bearing. 


Phe load capacity is given by 


2 . r^e 


W = vp^ rf + / 2iT rP dr 


which simplifies to the following equation after using equation (3.40) 


W = 6T1Q / 


(3.42) 


r. IT 
1 


Similarly the equation for frictional torque can be written as 


T = 2'n-n^ / 


r 3 


(3*43) 


where S3 is the rotating angular velocity. 

ITow, for a rough hyd 3 rostatic bearing, the fzlm—thicJmess is assumed 
to be approximated by 
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M 

I 

d=i 


H = h [ 1 + j ^ QQQ ^ j 


a. 

2 


where h is the nominal fiha-thidcness, e. are the amplitudes 

<3 

are the wave length of the rou^ness waves. 


(3.44) 


and 2tt a. 

0 


In the case of hydrodynamic lubrication and the 

h ’ 

function can be expanded and the various expressions mentioned before 
can be approximated as follows : 


griO 1 U J COS X. t 

P, fln-l - 3 t ^ / a_ 


dt 


2 

M e 


-A* . 1 cos 2X. t 

+ 3 I -|{ln-i+ / — dt }] 

j=1 h ^ X ^ 


(3.45) 


6nQ r 2 M e . 1 

W = 5 5“^ [■■■5 3 1 / t cos X. t dt 

h^ 2 3=1 h X 3 


M e 1 ,2 1 

3 I -4 + / t cos 2 X. t dt }] 

j=1 h^ ^ X ^ 


(3.46) 


T = 


pTTTiJJ r^ . ,4 M 1 •! 

' - I I cos X^ t dt 


h 


[ 


0=1 


2 

e. 


o' I “^ + / cos 2 X. t dt}] 


j=1 h 


r. r 

1 e 

where r = t, X= “ , = ~ 

8 0 


(3.47) 


Star fairly rough surfaces, the wave lengths of the rou^ness 
waves are sxifficientjy small and a. tends to zero. In such a case 

tends to infinity and all the integrals [see appendix III] in the 

<3 , , 

prevxoiis expressions sre jssero* Hence the expressions for detiexoiiniii^ 
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■fclie flow fLiiXy load capacitiy and frictional iorq^ue ars given respectively 
as follows ; 


6tiQ - 1 

— ^ In - {1 + 3 
TTh^ ^ 

I -i) 

3=1 h^ 


(3.48) 

3UQ / 


a 

2 

s. 


e 

^3 

(i - x^) 

[1 + 3 I 

-4 ] 

(3.49) 

h 


3=1 

li 


irnf^r^ 


1 ^ 

*1 1 

s2 


e 

2h 

(1 -x^) 

_i 1 

2 ■* 

A 

(3.50) 


It can be remarked here that, as the amplitudes of the 
rou^ness waves (roughness parameter) incr§§ee ,, the load capacily 
and dCfictional torque increase for a given flow flux. However, the 
flow flux decreases as o increases for a given inlet pressure. 

APPMDIX - III ; EVALUATION OF IXfEGRiiLS 


It can be shown that the following integrals tends to zero as 
X. tends to infinity 
1 cos X. t 

(i) / — 


sin X. sin XX. 1 cos X- t ,, 

— 1 1 + / — 

^ X b ^ 


(3.51) 


Since sin and cos ^1* between —1 and +1 for all values of X^ 

and are therefore bounded and hence 


sin 


X, cos XX. t 

— ^ i. ->-0, — i-— ^ 0, as X. ^ » f 

b b ' 

and the above integral tends to zero as X. tends to » . 



(ii) Since 


1 


/ t cos IV. t dt = 

X ^ \ 


sin X. cos X. sin X. cos x. 




-i - X 


JL . 






JL 


(3.52) 


and therefore by above argument this tends to zero as X tends to « . 

D 

1 3 

(iii) Similarly / t cos X. t dt also tends to zero as X. tends 


to <» . 


APEEKDIX-IV 5 DETERMINAKOI OP Ajs AUD B?s 
^ I3 

Hie coefficients A= ’s and B. ’s in the equation (3.4-a) and 

J. J 

(■Bi.4-b) are so chosen that the mean square errors 

■ tt M id V ^ i 2 

I ihii] } and E{[Q I b r] are minimum, 

1=0 3=0 

the conditions for which are as follows : 

M 

E [(Q^ ) - ? I A, ECif"'^)] =0, (v = 0,1,...M) (3.53-a) 

^ ^ i=0 

E [(Q 11^) f B eCh^"^)] = 0, (y = 0,1,. ..,B) (5.53-b) 

y 3=0 ^ 

The equation (3.53-a) gives a system of (M+1 ) equations with (m+ 1 ) 
unknown A|s piovided ^(Q^^h'') and e(H^^^) are known. 'Bir considering 
the density distribution function of h^ to be noimal with variance o, 
one can determine the values of e(iF) for any m [ see appendix 1] . 
Hie values of E(Q H* ) can be determined only when the correlation 
functions between and H are known for all v . 
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Since these correlation functions are not known we can 

deteimine them approyinately for a given bearing system as follows. 

For a given bearing surface, the roughness profile can be traced 

by a Talysurf and thus, for each point of the bearing surface, a 

corresponding film-thickness H can be determined. Ihen some 

representative values of H say are chosen. Ihe 

■^0 

fluxes jQg, . • . jQjj. corresponding to each of these’ representative 

0 

values of film-thickness can be determined from the cases of bearings 
with smooth surfaces having film-thicknesses 

0 

assuming tliat a given rough bearing is equivalent to number of tiry 
bearings with smooth surfaces joined together, we can write 
approximately as follows : 

S(Q H^)=,|- f Q^. , U = 0,1,..-,M) (5-54) 

i=1 i 

Similarly, the values E(Qyi?' ) for (y = 0,1,...,e) can also be found. 
Thus, when e(h''), E(q^h'') and E(QyH^) are knovn, all the coefficients 
Ajs and Bjs can be deteimined approximately « 



LIST OF HOTATIOIJS 


b 

E( ) 
E 




8 

kl 


L 


M,N,ia,n 




Q,Q. 




width of the bearing 

expected value or statistical mean 

frictional force 

nominal film thickness 

stochastic film thickness 

total film thickness 

step position 

Length of the step bearing 

real numbers 

pressure 

flow flux 

radial coordinate 

inlet radious 


T 

U,V 

W 


2:,y,z 

a 

h 

2 

a 

y 


torgue 

surface velocity components in x and y-direction 
load capacity 
coordinate system 
suffix 

lubricant viscosity (constant) 
variance 

friction coefficient 


K 


random variable 
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CHAPTER - IV 


EPEECTS OP VISCOSITY VARIATION IF STOCHASTIC lABRICATIOF 
4.1 IFTROnJCTIOF 

In the previous Chapter, a stochastic form of Eeynolds 
©(juation applicable to rou^ bearings has been derived bj"" considering 
the flow of an incompressible lubricant, having a constant viscosity* 
As described in Chapter — II, viscosily can vaiy with temperature and 
pressure etc., it is desirable, therefore, to consider the variation 
of lubricant viscosity across as well as along the fluid film while 
studying the bdhaviour of lubricated systems. One of the approaches 
to study the effects of viscosily variation was suggested by Tipei 
[1962 ]by considering a typical viscosity film-thickness relation. 
Tipei and Fica [1967 ] , Tipei and Degueurce [1974] have studied the 
effects of viscosity variation in the cases of thermohydrodynamie 
problems and shown that the solutions obtained for vLscosity-film 
thickness relation, considered for plane or curved surfaces with 
convergent film, are quite satisfactoiy, Perfomance of viscosity 
variation across the fluid film with various profiles have also been 
studied by Quale and Wiltshire [1972] . It was pointed out that the 
coefficient of friction reduced significantly if the viscosity of the 
lubricant is made to vary across the thickness of the fxlm. 
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Keeping in view of these, in this Chapter, a generalized foim 
of Reynolds equation for stochastic lubrication is derived by 
considering viscosity-film thickness relation as suggested by Eipei 

[1962 ] , In particular, the case of one-dimensional rough slider 
bearing has been studied. 


i\irther, a generalized form of Reynolds equation for stochastic 
lubrication for multiple layers of lubricant having different viscosities 
between two rough surfaces is derived by assuming Ihe corresponding 
Reynolds equation as derived in Chapter II, without slip, Ihe corres- 
ponding cases for longitudinal and transverse, one-dimensional rou^ness 
are given. Ihe case of rough hydrostatic bearing with three layers of 
lubricant with different viscosities has been studied as a particular 
case. 


4,2 viscosimr yariaiior adckg ihe rih 

Consider the flow of an incompressible lubricant having 
variable viscosity between a stationary rou^ surface and a moving 
smooth surface [see Pig.(4.l)-]. following Christensen [1969-70] it 
is assmed that the equation governing the pressure in the thin 


film is given by 


3. . £. (£ i|) = 6U P + 61 ^^ 


(ii_ , , 

3x n 3y 'n 3y' 


9z 


ay 


(4.1) 


The following viscosity film liiickaess relation is considered, 
Tipei [1962 ], Tipei and Degueurce [1974] , 
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where the value of 'q' depends on the operating conditions of the 

H 

bearing. Since ^ _ 1 , it can be noted that the viscosily decreases 

1 

along the film- thickness and it is hipest when q == 0 and lowest when 
q = 1. 


It is observed that the pressure in the fluid film is a stochastic 
parocess as the film-thickness function H is a stochastic process and 
can be determinod by averaging process, Papoulis [1965] and P^lenik 
[1968] . Proceeding as in the previous Chapter - III and using equation 
(4.1) and (4.2), we have 


9x 




E [hr H- 


■3-q 3P 


9x 


] 


ay 


E [h^ E 


q „3-q aP 


ay 


] 


6n u e(h) + 6n 

0 9x ^ c 




(4.3) 


where, 

<?o 

e(s) = J .s e(s) ds. - • • - . ■ - . 

— CD 

and P(s) is the probability density distribution function of the 
stochastic variable s. Ihe film-thickness function H can be written 
as ' 

H = h(x,y ) + hg (x)y, 5) (4.4) 

as described in Chapter - III. 

As before, it is assumed that the flow fluxes vdiich are 
stochastic process can be written in powers of H, Sveshnikov [1966] 


1 
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Lin [1967] , Christensen, et. al. [1975] 

.-n TT 


Q -••5 H 

X 2 


91 _ 5 H ^ ”1 


I2n dx 2 




3P _ U 
9x ~ 2 


M 


I A (4.5-a) 
i=0 


3 ^Q. p.3-q 

n - 1 H „ ^1 ® 3P Y ? i 

- O h - -7;:^ ^ _ H _ ^ - X Y B 

'2^0 S7 2 -io 


(4.5-b) 


y 2 " I2n 8y 2 i;in_ Sy 

where s and s are non-random functions of x,y,U and V or constants 
depending upon the type of roughness, roughness slopes, Tni pimiTn} film- 
thickness etc* and may be deteimined theoretically or experimentally 
[see Chapter III ] . 


Poliowing the same procedure for an analytical evaluation as 
discussed in the Chapter - III and using equations (4.3) , (4.4) and 
(4.5)1 equation governing the fluid pressure can be obtained as 


follows : 


E(lf ) h^ S(h“) h^ 

_[__^™e(p)] [-__)^E(P)] 

6 n. u [ I a^e(h^) + — tE( 4 T) - I -a,.e(h^‘" 5 ’^)>] 


ax 


i=0,i/4n 


g(j^-3+q) . jj2-q 


+ 6n V|r [ f B.E(I# ) + — {E(-^)- I . B e(h3"^^)}] 

0 ay e(h^^^) i=o,3;^n3 


,i=0,3i4n 


(4.6) 


Equation (4.6 ) represents the generalized form of R^nolds equation 
for stochastic lubrication applicable to rough bearings with viscosity 
variation. When q = 0, this equation (4.6) reduces to the same form 
as discussed in Chapter - III [see equation (3«8)] •« 



In the following, the case of an infinite slider bearing with 
transverse, one-dimensional roughness is considered. 

4.3 imNIIE SLIDER BEifflING 


Consider the case of an infinite slider bearing with transverse, 

one -dimensional roughness as shown in ilg. ( 4 . 2 ). ng A. = 0 for 

ij4i 

all i and m = 0 in the generalized equation (4.6), the E^nolds equation 
applicable to this case can be written as : 


h 


<-i 5 ) 




1 




E 


Integrating equation ( 4 . 7 ) with respect to x, we have 

is [E(^)-0E(^)1 

1 

where C is an integration constant. 


(4.7) 


(4.8) 


Integrating equation (4.8) further and using boundary conditions 

E(p) = 0 at X = 0 , 

E(p) = 0 at X = I. 

we get, 

6ti^ U X . 

0 r r » 


e(p) 


J [E(-|::j) - 0E(^)1 ax 


h^ 0 E 
1 




where 


e(-^) dx 


0 E 


(4.9) 


(4.10) 


(4.11) 
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•me expected value of load capacity is obtained by integrating the 
eiprosaion of E(p) over the bearing surface area as follows ; 


E(W) = / sWaotw-/"' X I-Wp)] dx 

0 


0 


( 4 . 12 ) 


This gives, 


6ri U I 

e(w) = - ■ ° 


_/ X te(^) - 0E(^)] dx 


(4-13) 


The expected volue of shear stress on the runner is 


= [ 4E(^) - 30E(^) ) (4.14) 


The friction force over runner is obtained by integrating the above 
expression of E(t ) as follows : 

E(p) / [4E(-^) - 3CE(~^) ] dx (4.15) 

h'l 0 

.%■ using approximate values of various quantities involved 
in the above equations (see Appendix), the expressions for E(w), E(p) 
and E( v) can be calculated apprcximately as follows for 0 q 1. 


Por 0 < q < 1, 


u 

0 


rr 1 . 1 1 , 1 1 t 

q (i-q ) g q ( 1 -q ) 

1 1 

*" ^q ^ (i-^)( 2 -q) rq (' 2 '-^') C2 (i-q) c ^ ^ 


, -ilTii +i-} 

^ h^ 

_ a^ii^ .11^+^ }] (|-)2+a^(0) (4.16) 

’ n hi by a 

1 1 1 
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■gjr E^l^m [1 , _I_ _,Lo (I. , J_',, 
71 D ,-q ii-q) — )] 

o K ' S hi 


[2(l-q) (^-^) -|o^ ( 2 -q) (i-Z?)] (f)^ +^"(0) 


h^ e; 


E^ h, 


where, 


C 


(2-1) t * “'‘(0) 


lij h, 
1 1 




fff h^ ^ h^ 

For <1 = 0, 

^[wj = = [m h^ - (i - ■^) - ~ (1 - I 


(4.17) 


(4.18) 


t(i-f 


h^ h:j 


,-3 ^'■‘ "h 

hi hi 0 


4 -3Ci (1 -^) 

n 11 r 


(4.19) 


+ 2 (l - -^ 5 ) - 3 C^(l - ^) (^)^ + a''( 0 ) 


h 


where, 


(1 _-^) + (1 {^f + o^{0) 


h 


= 2 
0 


For <1 = 1, 


1 


1^5 K 
hi 0 


(1 -■“) + 3(1 -^)(|-)^ +'"^^0) 

h, o 


h 


1 


Mil 


‘1 


^,n h 


Jin h. 


m 


6n,u 


<10 XX * .4 — .r-Tf ^ 

^ (1 _ 1. _ _^) + Ci(r r~l] 

h. h, hi hi 


,.a.t(r-i.)2-Ci(i -^-%)3(r) - 
2h, h, -ill I 1 -, ° 


■ #(0) 


( 4 . 20 ) 


( 4 . 21 ) 


( 4 , 22 > 



^= 4 ^ = ( 4(1 

° ' ' ‘^1 


3C 


1 /I 




2 'h S5' '!i‘ 

“l *'1 ° 


where, 

iin h 


1 +1 ri. 1_\ (O s2 

~ 2 'h ^ + cr (o) 

N N K ° 

1- -L.'\ ^ \(^ ^2 4/ s 

h^ h^ h^ 0 


C, 


(• 


(4.23) 


(4.24) 


!]}he oocifxcxexiu ‘j.l Xric*fcxon. on "the ixinner cen be d.e£Liied. xn ijhe non* 
dimensional form as follows : 



(4.25) 


iTo see the effects of roughness and Tiscosiliy variation the 
expressions for ]3(w), E(f) and E(y) are plotted in the Eigs. (4.3), 
(4.4) and (4.5) for q = 0.0 , 0.5 and 1.0 respectively for different 
values of . It can be seen from Eigs. (4.3), (4.4) and (4.5), that 


0 

the load capacity and friction force increase for all values of q 

0 ^ 
as ^ increases but the coefficient of friction decreases as — 

0 o 

increases, further, it can be seen from Figs. (4.3) and (4.4) that 

the case corresponding to q = 0.0 attains higher load capacity and 


more friction force than that of the cases corresponding to q = 0.5 


and q = 1.0. However, from Eig^? (4.5) it has been seen that the 

coefficient of friction in the case q = 0.0 increases up to a critical 

value of ~ ^ 0.18 and then decreases as — increases from this 
o 0 

critical value. 












4.4 VISGOSITT YARIATIOII ACROSS !EHE PITM 


In the previous section 4.2, a generalized form of 
Hsynolds ecjuation applicable to rough bearings is derived by 
considering the variation of fluid viscosity along the film. In 
this section, we derive the stochastic Reynolds equation vdien 
viscosity varies across the film as a step function. As derived 
in Chapter - II the Reynolds equation in the two-dimensional form, 
without slip applicable to three layers of fluid can be rewritten 
as follows, [see equation (2.40) and Rig. (4.6)] 


9x 9y 9y'' 2 9x at ^ ’ 


whe3re 


h-’ , (h+H)^ - h" 

^'1 — 


G 


(h+H) 


Fluxes in x and y direction are given by 

R. 


Q, 


X 


u 

2 

R. 


G 


Si 

12 9x 




12 9y 


(4.26) 


(4.27). 


(4.28) 


Considering the bearing configuration as shown in Rig.. (4.6), the 
film-thickness in this case can be written as follows • 


H = 1^0 


h+H = h + h„ + hg = h^ + h 


(4.29) 
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where h^ is the stochastic part of the thickness ana ^ is the 

noiminal film ■ '-iLdcness for upper and lower layers of lubricant 
having viscosity n^. 

Taking the expected value of equation (4.26), we have 


f- E[f^ , 
3x 1 ax 


O [J, ai] ns 


ay 


1 ay “ 2 + at + ^g) ^ (^*50) 


To ©valuate the various terms in the above equation in the case 

h 

oi hydrodynamic lubrication, the ratio may be assumed to be 


small, i.e. then have, 


n 


S<hn + hU = h^ E(1 + = h_^ i h_^ E(^) 


(4.31 ) 


n 

Assuming that the probability density distribution for h to be 

s 

normal and having variance a , the expected value of E(h + h ) can 

n s 


be written for r* « 1 as follows ; 


E(h + h ) = h 
^ n s n 


(4.32) 

By considering one -dimensional roughness pattern when the roughness 
direction is parallel or perpendicular to the sliding direction, as 
suggested by Christensen [1969-70] and Christensen and Tender (19693 , 
the Reynolds equation in the following cases are derived from equation 
(4.30). 


4.5 l ONOiroBIIlAl. ORE-DIMMSIOR.AI RCTJGHgESS 

In this case, roughness is assumed to be ef the form 
of long , naxT’O" ridges and valleys running in the direction of 
sliding and the film thickness (h-Rl) is represented by, 
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h 

where h is the stochastic part of the thickness and — is the 
^ 2 

nominal film ■■■ bickness for upper and lower layers of lubricant 
having viscosity . 


faking the expected value of equation (4.26), we have 


L- 

dx 


If ' If ® ^1 w ^ ^ . hj ] (4.30) 


2 32: 


9t 


To evaluate the various terns in the above equation in the case 

h 

of hydaxidynamic lubrication, the ratio — may be assumed to be 

o n 

small, i.e. « 1, we then have, 




(4.31 ) 


n 


n 


Assiiming that the probability density distribution for h to be 

s 

normal and haTi:ig variance a j the expected value of + iig) can 
be written for ^ « 1 as follows : 


E(h + h ) = h 
^ n n 


(4.32) 

By considering one— dimensional roughness pattern when the rou^ness 
direction is parallel or perpendicular to the sliding direction, as 
suggested by Christensen [ 1969 - 70 ] and Christensen and fonder [19693 , 
the Reynolds equation in the following cases are derived from equation 
(4.30). 


4.5 lOHSrroBUfAl. ORE-PBiEH'SIORin:, ROKaPTESS 


In this case, roughness is assumed to be of the form 
of long , narrcvy- ridges and valleys running in the direction of 
sliding and the film thickness (h-ffl) is represented by, 
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h+H = hj^(x,y,t) + h^(y, g) 


( 4 - 33 ) 


where h„ is nominal film thickness and h is the stochastic part 
of the film- thickness . 

itollowing Christensen [1969-70 ] and assuming that ~ is 

9x 

gp 

a variable with zero variance, and can be considered to be 

3x 1 

(approximately) stochastically independent quantities. Shen we 
can write. 


[E(P)1 


(4.34) 


Also as flow flux in the y-direction 


dn nn'h vpi tv 


can be found as follows ; 

Considering 

I’^|| = qy(x,y) ( 4 . 35 ) 

which do not vaiy randomly, we write 

and, 

n _ 3(p -^) = __1 2 >_ g(p) ( 4.37 ) 

Putting the expressions (4.32), (4.34) and (4-37) into the equation 
( 4 . 30 ), we get 

fe (Kf, ) ^ E(P)] . ^ 1 ' f ^ It < V 

(4.38) 
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This is the Eeynolds-iype equation for a bearing with three layers and 
applicable for longitudinal, one-dimensional rou^ness, 

4.6 IKMSVERSE, QgE-hlMEffSIOErAl ROUGHUBSS 

In this case surface rouglriess is assumed to have the 
form, of long , narrow ridges and deep valleys running in the y-direetion 
(perpendicular to the sliding direction). Then, the film-thickness is 
given by 

h + H = h^(x,y,t) + hg(z, C) (4.39) 

Following Christensen [1969-70] and carrying out a similar calculation 
as in the previous case, we get 


fe ' ETiTiTT « ■"ly iF 

(4,40) 


eCg^f^ 


which is a Eeynolds-type equation for a bearing with three layers of 
lubricant and applicable for transverse, one-dimensional rou^ness. 


Similarly the stochastic equations governing the pressure 
in oases of other forms of rou^ness can be derived, Christensen [1969— 7*^ 

4.7 HYDROSTATIC BEARCTC \?IgI TliREE LAYERS 

Consider the flow of three layers of lubricant in a 
rough bydrostatic bearing as shown in Pig. (4*7). The Reynolds 
equation for longitudinal roughness in this case can be written as ; 

|j:trE(Fp|jE(y),]=o (4.41) 

where P.j is given by equation (4.27), 
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Integrating equation (4.41 ) and using the deiinition of flow flux, 
we get, 

E(q) = -%TrVE(i'^)|^E(T) (4.42) 

The ejqpressio ns for the pressure and flow flux can be obtained, by 
integrating equation (4.42) and using Ihe usual boundaiy conditions 
E(1>) = E(p. ) at r = r. , and E(p) = 0 at r = r , as follows : 

”i.TTE(p(j (4.4J) 

The Expected value of the load capacity can be calculated as usual 
and we get, 


e(w) = 



(4.45) 


Considering equation (4.27), for hydrodynsaiic lubrication as considered 
before, E(E^) can be written as ; • 




t? E 

12r)2 12n^ 


<[1 


12n, 


+ ■ 


3(f)' 1 

n 


- h" 


(4*46 ) 


I2n, 


It can be seen from equation (4.46) that e(e^) decreases 

as b. and n_ increase and E(e ) increases as (^) increases. Hence 

■^n 

for a constant flow flux, the expect ei value of load capacity increases 

as n and n increase and E(\y) decreases as ^ increases. The last 
' . ‘ - 'n 

result is the sane as discussed in Chapter III for this case ishen m = 3 
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4-. 8 RESUITS MD PISOISSIOMS 

In this Chapter, a stochastic E^nolds equation applicable 
to rough slider bearing is derived by considering the variation of 
fluid viscosity along the fLlm- through viscosity-film thickness relation 
and a particular case of rough infinite slider bearing is discijssed. 

It has been shown that the load capacaly and friction force decrease as 
the index parameter ^q*” increases from zero to one while the load 
capacity and friction force increase with the rou^ness. It has also be 
seen that the coefficient of friction decreases as q increases upto 

a critical value of ^ * 0,18, and after that it increases as q increases. 

0 

Further, the stochastic equations applicable to the 
longitudinal and transverse, one-dimensional roughness are also deidved 
by considering three different layers of lubricant with different 
viscosities. Is a particular case, a rou^ hydrostatic bearing is 
considered. It is shown that the load capacity increases as the 


viscosities of lubricant layers increase 
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-APPMDIX 


In the case of hydrodynaEic lubiication, the ratio ~ may 

ll 

be assumed to be small where ^ «1 and one can write, 


h“ = . 1,“ [1 + B + -alfii (^f * ...J 


(4.47) 


Assximing that tbe probabulity density distribution for the stochastic 
variable h^ to be normal and having variance a , the eapected value 
of eqmtion (4.47) can be written for ^ « 1 as follows ; 

B(H°) * h^ [1 + -SlaiLl (|)2 + a4(o) ] (4.48) 

Similarly, one can write 


Por a slider bearing as shown in Hg. (4.2), 


(4.49) 


h = h^ - HOC 


(4.50) 


where m-the angle of inclination between two surfaces. We can have, 
dh = -m dx ( 4 . 51 ) 


Now, using equation (4.5l), we can write 

. X 

dx = 

m 


h 

/ E(h'') dx = i J ^ E(B^) dh 

0 h 


and ll h 

J E(lf) dx =;i/ ^ E(H^) dh 


(4.52) 

(4.53) 


0 


In a similar manner, other expression can also be written fer various 
values of n and integration can be carried out in the same we^ as 


discussed in Chapter III. 



LIST OP FOTATIOirS 


E( ) 


P 




H 

L 


P 




r,e,z 


T 

U,V 

W 


3!:,y,z 


n 




K 


expected value or statistical meac. 

friction force 

nominal film thickness 

stochastic film thickness 

total film thickness 

length of bearing 

real numbers 

pressure in the film of the lubricant 
inlet pressure 
flow flux 
povrer index 

cylindrical polar coordinate 
inlet radius 
outer radius 
torque 

surface velocily components in x and y direction* 

load capacity 

rectangular coordinate 

variable viscosity 

viscosity of the base oil 

random variable 


a 


varience 
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CHAPTEE - Y 


A riEW THEOIff- OP MIXIID HJBRICAHOII 
5,1 IFTROCTTCPIOU 

In thki previous tvvo chapters, generalized foims of Reynolds 
aquation for hydrodynamic lubrication applicable to rou^ bearings 
have been derived. In the case of rough surfaces the hydrodynamic 
conditions are maintained by using a high viscosity lubricant such 
that the film-thickness is greater than the heights of roughness 
asperities. If it were not the case, asperities penetration may 
occur and boundaiy lubrication can exist under extreme condition 
where bearing surfaces are in continuous contact at the tips of the 
asperities which support the load. Hae formation of oxide film and 
adsorbtion of the lubricant on the surfaces play an essential role 
in the lubrication mechanism (such as reduction of friction) but the 
viscosily of the lubricant which m^ be present in between the groves 
in a minute quantity does not take important part in the process, 
Bowden and Tabor [ 1950, 1964] . 

When the fluid film, which seperates the two surfaces in 
a lubricated system becomes thin enough such that the surface asperitie 
begin to interfere, a situation between hydrodynaaic and boundary 
lubrications, known as mixed lubrication arises. In this case, load 
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is partly supported by fluid film and partly by asperity-oontacts. 

The frictional force depends upon the interaction of the surface 
asperities with the fluid film and ihe total friction would arise 
partly due to boundary lubrication conditions associated with 
asperity-contacts and partly from bydrodynamic friction, JUller [1954' ]. 

Lenning [i960 ] has pointed out that the 'transition from 
bovindary to mixed lubrication occurs by increasing base oil viscosity, 
decreasing surface roughness and decreasing applied load. Dobny [1964] 
has also considered a mathonatical model to study the velocity of 
transition from boundary to mixed lubrication and showed that the 
velocity of transition depends upon the surface rou^ness, lubricant 
viscosity, load and external pressure. A model for the contact of 
hydrodynamically lubricated machine eomponents was also developed 
by Eiomson and William [1972’ ] from which asperity to asperity contact 
load could be obtained. 

Christensen [1969] and Fowels [1971] have considered the 
asperiV deformation together with hydrodynamic lubrication effect 
for elastohydrodynamically lubricated bearings. Christensen [1972 ] 
has also presented a mathematical model for mixed lubrication tdrich 
arise from an interaction of the surface asperities with hydrodynamic 
fluid film between the sliding surfaces. It has been shown that in 
mixed lubrication regime, friction is mainly controlled by lubrication 
properties of liquid-solid interface where as load is entirely 
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cdntrolled by the hydrodynamic properties of the bearingi fsao and 
Tong [1975] presented a model for mixed lubrication by representing 
the surface asperities by short plasto-elastic cylinders with 
spherical tops where the heights of asperities afe assumed to have 
a Gaussion distribution. Solution of Reynolds equation was combined 
with boimdaiy lubrication theory or the micro-elastohydrodynamio 
theory to study the condition of mixed lubrication. 

In. this chapter, a new theory of mixed lubrication is 
presented by considering the surface asperities to be randomly 
distributed fins or cones which penetrate through liie fluid film 
and partly may be in contact with the asperities of the opposite 
surface, IQiese asperities are considered to form a net work of 
interacting zone corresponding to each surface, through which the 
lubricant flows as happens in the case of flow throu^ channels with 
singular resistances or some what as in the case of a flow throu^ 
porous matrix. Such situation m^y also arise when the additive 
present in the lubricant, such as long chain polar compounds (fatty 
acids), get adsorbed on the sliding surfaces with their free tail 
penetrating into the film, and thus forming a thin interacting zone, 
Askv/ith, et. al. [ 1966 ] . 

Ror the sake of generalily, it is considered that these 
interacting net v/ork zones are seperated by a thin lubricant film 
whose viscosity may be different from that of the viscosities of the 
lubricants present in the two interacting zones due to surfactants, 
etc., Askwith, et. al. [1966] , Rein and Kreuz [1966] , Ihus, the 
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effective lubricated region is divided into two interacting surface 
net woite zones and a middle zone through which the simultaneous layers 
of fluid with different viscosities flow. In these interacting net work 
zones, modified forms of Havier-Stoke’s equations are assumed and 
which can also be thought of as generalized forms of Darcy's law 
applicable to porous matrix, Brinkman [1949] . However, in the 
asperating film, the usual hydrodynamic equations are considered- Keeping 
these points in view, a generalized form of Eeynolds equation for mixed 
lubrication is derived and particular forms applicable to various cases 
are deduced. 

A typical situation where this theoiy can be applied is the 
case of synovial joints under impact loading corditions. Ihese joints 
consist of two poroelastic surfaces called 'articular cartilage', with 
a fliiid lubricant film between them. Ihis fluid, known as synovial 
fluid, contains hyaluronic acid molecules, which changes its viscosity 
and also acts as a boundary lubricant at cartilage surfaces under 
extreme conditions of ■ operation, Dowson [1967 ], Dowson . ■ . 

et, al. [1970] ,rGhandra [1-9753 jShukla aiid Chandra [1975 1 > ■ 

Ihe superficial tangential zones of the Cartilage and the region of 
asperities peneration in the synovial fluid film may be considered 
as the interacting surface net work zones while studying ^novial 
joints by the help of the theoiy presented in this chapter, IiJittal 
and Millington [ 1971 ] , Mow and Eedler [ 1974 ] . 
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5.2 BASIC BSUAIIOIS 

Consider the flow of three incompressible layers of 
lubricants in the two interacting 2 Dnes and the middle separating 
layer which mi^t exist atleast in the region of hydrodynamic 
lubrication. Bor example, in the slider bearing, mixed lubrication 
can exist at the trailing edge while at the leading edge hydrodynamic 
lubrication can preTail, Christensen [1972] . When the mixed 
lubrication condition exists eveiywhere in the bearing contact 
area, the flovir in the separating middle layer can be ignored. In 
general the physical situation and the coordinate systems described 
here are shown in jBig. (5.1 ). 

using the usual assmptions of lubrication theory and 
considering the same pressure gradient in the three layers, a genera- 
lized form of Reynolds equation is derived as follows ; 


In the interacting region I. ^ z £ + H£), fluid 

motion is considered to be governed by the following modified 


equationg of motion, Brihkmann [1949] 
2 


9 u^ 

9P 

Til 

^1 . 2 
9z 

“ ^ 

“^x 

,2 

9 V. 

9P 

h. 

'’l ,"2 
9z 

“ iy 

<t> 

j 


(5. 1-a) 
(5.1-b) 


where ~ and ~ measure the degree of resistance to "the flow and can 
^x '^y 

depend upon the characteristics of the net works which in tern may 
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d epend upon the type of asperity, its distribution, heists, shape, 
molecular structure, crystalograply of adsoibed molecules etc. In 
the case of a synovial joint, these depend upon the characteristics 
of the superficial tangential zone of the cartilage surface and 

1 

the molecular structure of the hyaluronic acid. In fact, when — > 

-i- increase the residence to the flow increases and they are zeio in 
4iy 

the case of no resistance. 

In the region II (h^ + ^ ^ + h), the fluid flow 

is governed by the usual equations of motion 


3 ^2 3P 
^2 ~~r “ ai = ° 

dz 


(5.2-a) 


3 v 3p 


(5.2-b) 


The thickness of this hydrodynamic film layer be negligible in 
that part of the bearing vhere mixed lubrication exists. 

In the region III (H. + H +h<z<Hp), similar forms of 
equations of motion as in region I are considered to study the flow 
characteristics in hae second interacting zone. 


3 u^ 3P n3 


5 3z^ 


3x ” 4 


(5.3-a) 


3 v„ dP 


(5.5-b) 


The equation of continuity can be written as follows : 
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M + M +M = 0 

3x ay az 


(5-4) 


Now for the region I, integrating equations (5.1-a) and (S-l-b) with 
respect to z and using the following boundaiy conditions 


u^ = at z = Ilj 


^1 “ ’ ^1 ~ ^ ~ ^ 


the expressions for fluid velocities can be written as 
z-H^ _ 2 - (h^ + %) 


Ui =U. 


I z - (H + H ) 

J. _ sinh ] L. 


sinh^ 

— sinh 

, . T. 

T-H-. sinh — 

/^y 

z - H z - (H^ + H^) 

sinh sinh 

ch - ^ /T 

.li ri :i 1 

n I * TT. 


(5.5) 


■1 f (5.6-a) 


(5.6-b) 


In the region II, integration of equation (5.2-a) and (5-2-b) with the 


botindaiy conditions 
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Ug = Uj, , Tg = at z = + % 

( 5 . 7 ) 

Ug = Vg = at z = + h = Hg - 

give, 

-I p 9^^ 

“2 (H,« .h)] - 

u.-u. 

+ (-^) z + ^ [u^(H^tB^+h) - IL(h^+H ) ] (5.8-a) 

and 

12 

^ U - + (H,« ) (H,fH rt)]- 

7 . -- 7 . 

+ (-^;~) z + ^ [T^(H^+H^4ii) - 7^(H^+H^)] (5.8-b) 

For the region III, again integrating equations (5.3-a) and (5-3-b) 
and using the boundaiy conditions 



sinh 



siiih. 


sijQh 


^ ^ 


z - (H^ - %) 


V, a V. 


sinh r 


H^-z z - (H - H ) 

sinh — 4- sinh — 


[. 

^3 


1 J ay 


(5,10— b) 


To find the expressions of U^, V^, U. and 7., shear 




gas at 


interface i at z = H| ^ interface j at zi = 


= H + 


jj + b = 


ai*e ©juated as follow ; 


3^^ 3^2 

3z ^2 3 z 


3 ^^ 9^2 

3z "" ^2 9z 


at z = 




^2 9z ^3 3z 


3V 3V 

’^2 ^ ’^3 ” 


at z — 


•fiire® 


After using the expressions for velocity gradients 


regions in equation (5.11 ), we have 


^i^l " ^3 h “ ^1^3 “ ®1 9 x 


bo 3P 

-Ui ^ f2 = '^2 ^4 " ^2 ^ 


(5.12-a) 


(5.12-b) 
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Up ap 

7 f I «. V -»£ - V f » ~ — 

'^5 h ^ ^3 ®l 3y 


\ ap 

>Vi h = ■'^2 ^4 " ^2 9y 


where 


ho h^ 


fi 


f • = 


^3 “ 



tanh — ~* 

/T 

^2 ^ ^ 

1 

^ /T” 

h 

tanh 

Ay 

^2 ^!! 3 _ 

1 


tanh 

/T 

X 

\ ’^3 

1 

^/T ' 

tanh 

n^ 1 




sinh 


1 

^3 

% 

sihh 

1 


sihh. — — 

■'■S'l 


1 

. , 

sinli ~— 


(5.12-c) 


(5.12-d) 


(5.13) 
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g s ^ + 1 /^ tanh 
(^1 2 ^ 


2 47 

X 


g. = I H- /r tarii -i— 

2 y z/r 

J 


2 /<^_ 


V I .. ..■ - 

gi = -S + tanh — === — 

«2 2 y ^ 

y 


Solving equations (5.12), the ea^iressions for U . , V. , U. and Y. can 

-L X 


be obtained as follows : 


n2 . 9P 


"2 - (h^)‘ 


h ^4^2^ " ^^1^2 h 


(5.14-a) 


\ = 


il ig ^-h ^ 


■ Ho ' 3P 

rfv f ’ f ' + — ^ f ’Y ) - (s'f * + ff’ ) 7 — 1 

'■^12 3 h 4 2^^ 2 h ^2'^ 3y ■* 


U. = 
D 


f f - (P^f 

^1^2 ^h ^ 


9 \ / ^^9 \ 

[(^2=1^4 + -g- -3"^^ ) - (§2^1 + ~ 3 


(5.14-b) 


^3 = 


fffj _ (— £)' 

^r2 ^h ^ 


h Po 9P , 

[ (Y f ’f ! + — f *Y ) - (fi-'f ’ + — eM — ] 

*■^ 2 14 h 3 V ^^2 1 h ^1^3y 


itirther, in. this case, the flow fluxes can be defined as follows : 

H H. 4ii, H -H H 

Q = /2ud3=/l Vdz+/^ dz + J*^ u,^ dz (5.15) 
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■V =/ 

*7 tr 




V dz 


H, 


1 


H -H H- 

f dz + f Vp dz + / T_ dz 

H, 


After using the exprsssions for velocity in differejot regions 
appTOpriately, these equations give, 

— % r— ^f 

^ ^ 2 ^ 2 »c ^ 


-5 


% 


- + — (Hff - 2 tanh 

12n2 ^ ^ 2 ^■ 


^TT r"^ 

+ — (h^ - 2 ^ tanh i-)] — 

^3 ^ ^ 2 /r ^ 


Qy = tanh 


% H 

+ Vp >C tanh — i— +Y g > + Y g ’ 


-l-^ + ^ (a - 2 /i" tanH -^) 

12 2 1l ^ 2 ^ 

y 

+ ■;*- (H„ - 2 v^ tarn — —)] — 

f y 2 


Integrating equation of continuity (5*4) with respect to z and 
the limits from z = to z = F^, one gets 


H, 


H. 


2 ^ 
3x 


Hp »v ®2 

dz + / ^ ™ dz + [w] = 0 

H 1 


(5.16) 


(5.15-a) 


(5.16-b) 

taking 


(5.17) 


Using, well known fomula 



l-r-j- 


a , ‘’“2 

^ f(x,y,2!) dz * J ^ f(x,y,z) dz + 

9 x (L H 


f(x,y,H2) 


-1 . f (x,y,H. ) 


( 5 . 18 ) 


equation (3.1?) can be rewritten as follows : 


A- J ^ U dz + ~ ^ T dz - Ug gx 


_Y ---1 + V + [w] ^ = 0 (5.1 

''2 ay 1 sy 

substituting the expressions of and in equation (5.19) 


(5.19) 


No w o n 


we have 

3_[{ J- H. ^ (H,-2 tad. + ^(%-2 tanh 3 

3X 12n2 ^1 2 3 2 


(H5-2 tahh 
av *• I2n„ n. ^ y 


3y '• 12312 ^i-i 

_ TT tanh 

- '^1 3 x '• 1 X 


) + (Hf-2 ^ ■tanh 7-“)> ^7^ 


2 3/^ "3 


2^ 


_i^L— 1 _ TJ -— [H_ - /$” tanh — ] 

2 /4.^ ^x 


, -rr ^ f F. + tanh 

^ ^1 3y '• 1 y 


^ - I _ ir — [H„ - tanh — ^ 

2 3^ y 


"Is ("i * "3 *''3 * “s' 

This is a generalized form of Reynolds equation governing the 
pressure under the condition of mixed lubrication. 


(5-20) 



155 


In th© following, particular cases are derived from the 
generalized equation (5.20) and their one dimensional forms are 
studied, 

i. SHDEH BEAB.IFff. 

ii. HTDBDSTATIC STEP BEAEEFff. 

iii. SQUEEZE EIIM BEARIir&. 

It is pointed out here that in the above cases, we have 
also slu-died the situation where the condition of mixed lubrication 
arise either in part of the bearing region or in whole of the bearing 
contact area. In the region where mixed lubrication exists, we 
have calculated the load capacity and the friction force only due 
to the hydrodynamic flow and the calculation of the components of 
these quantities which arise because of asperity contacts has not 
been included here. However, this can be easily done as suggested 
by Christensen [ 1972 ] and Tsao and long [ 1975 ] , 

5.-3 SLIDER BEARING UEDER IICXED lUBEECATIOU : 

Consider the case of a slider bearing where the lower 
plate is aaooth and moving with a constant velocity U. !Ehe upper 
plate forms the aforesaid interacting zone due to asperities through 
which the lubricant with viscosity flows. It may be considered 
that the region 0 £x £5.^^ may function under the condition of 
hydrodynamic lubrication and the region i il is under mixed 
lubrication, Chilstensen [1972] . 
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^ choosing, 


= 0 , Hg = h + H , 


= u, ITg = = Vg = 0 , [w]“ =0 

equation (5.20) reduces to the following form in the region 0 <x ^ L 


a 8 U a r 

i_ rj. — 1+ £_ TO 1 = u __ JJ J 

ax ^ X 9x J ay y 3y •' 2 ax x 


(5.21) 


whe3:« 


^3 ^ M - tanh M 

p + X, 

X 


. ^ aM + tanh M 
hir r X X i2 


tanh 2M 


nn r 


i? I? , 

y 12T12 4 n^ ^ 


[ 


2 ri. a M + n<-, tanh 2M 
1 X '2 X 


- tanh M 


(5,22-a) 


, tt 2 ctM + tarih M „ 

+ M_ Lf r 

T- !_ - iTf J I 


tanh 2^y 

2n^ aMy. + tt 2 ianh 21^^. 


a M + tanh M tanh 2M 

(j = [ ^ S 

^ ^'^x 2n^ oM^ + i^h 2M^ 


(5.22-b} 


(5.22-c) 


, „ ,-iL_ 

« " a ^ 


, M ^ 

a’C 


and M^, may be called ’interference parameters' and they increase 


as the resistance to the flow increases. 


ime equation (5-21 ) determines the pressure in liie region 
0 < X < il. where the viscosities in the layer 0 y ^ h and h£, y^ h+H 
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aro taken to be different due to the presence of surfactant in the 
lubricant. However, if this is not the case, - Hg can be easily 
taken. ■ 

In the region 1. x £ ,the equation (5.21 ) deteimines 
the pressure with h = 0 and the functions 1?^ etc. are modified 
according. 

iUo see the effects of the parameter of iiie interacting zone 
on the characteristics of the bearing in a simplified form, the case 
of an infinite step slider bearing is considered vriaere the step is 
located at x = In the region 0 £ x <. ii, h = h^, it is considered 

that the bearing functions under hydrodynamic condition with inter- 
ferance of the interacting net work zone while in the region ^ ^ t 

h = 0, it functions under mixed lubrication with some flow through 
interacting net work zone corresponding to this region. Ihe physical 
configuration is illustrated in Pig. (5-2). 

In this case, Reynolds equation (5.21 ) reduces to 


dPi 

dx 


const. 


i = 1 in region I 
i = 2 in region II 


(5.23) 


Since the pressure at the step p is continuoirs. 


dp dp 

p = £. = - ( ^ - «^ ) 

^s 1 dx 1 dx 


(5.24) 


Using the contimily of flow fluxes in the two regions and remaabering 
equation (5. 15 -a), we o an write 








aisH 


BEAr-iric m)V.^ ^s^xed lube icat ion 
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^ ® \ dx 


dp 

dx 

dp^ 


1 + I (h+ H G^) 


(5.25) 




where and G^ are 


defined in equations (5.22) and, 


F. 


h’ (- 


- tanb 


0 4n| 


m: 


(5.26) 


X 


tanh 


solving equations (5.24) and (5.25), the following expressions 
pressure gradiaat can be obtained. 


for 


dp^ G(i^- \) 


and 


dx 


to 


u n. 


(I - «.) + Hi 


(H- »i) I',, +Hi 


, 0 <_x < \ 


, 1. X 1 *- 


(5.27) 


( 5 . 28 ) 


Using the boundaiy conditions 
"n ^0 at X ^ *0 

p^ = P2 s-t X = )lj_ 

= 0 at X = 2. 


P2 

the expresslo 


(5.29) 


lou for the load capacity can he eyaluated aa follows : 


\¥ 


2 . 2 - 
J ^ p dx + / P2 


% 


S. hp^ 


0 


"i 



(5.30) 


(,.P7) ana (5.28), this gives, 

On vising o^^^atloas 

U ZS,. ( 1^- ^i) r ^ -i— ■ ] 


f B 


[ 


( ii- ^i) 


itollovwing the usual procedui'’ , 

can be obtained as follows i 


shear stress 


at the moving siarface 


dp 


2 


i (h + H&3,) >5^ H 


^ V 

2an^ + '^2 


for 0 5 X 5 


To=|^o‘^^ H ^tanh 2M^ 


*^^2 , for S-i 1 3C < 


■2 " 2 


jTow the force 


of friction at the moving 


(5.31 ) 

(5.32) 

surface can be obtained by 


integrating the 


above expressions over the 


interval 0 < x < a. 


i*e- 


? = J ^ 1 
0 


dX + / T. 


which gives, 


U a- Z^) - 0^)1 

T ( S,- !l . )P^ + S!-i^'o 

n u( a- i^j_) 

'tanh 2I''L 


+ g 

H 


tarih 2M 

a . 1 ^ — j 


(5.33) 


2M 


^ see the various effects 

•+,r and the force of friction, 
load capacity 


n2 

of interference parameter 

it is note! frcm the 


on the 



appendix that tho :function [h - H(g - G^) ] is positive and increases 

w A 

tanJa 2M 

03 m (a -iL*-u.) increases. i\irth©r, tie functions, f , JV, — rr? — ^ 

X 2 ^ X 0 

3C 

decrease as VL increases. 

A 

!Ehen, it can be seen from equations (5 30) and (5.33) that the 
load capaci-ty and the friction force increase as increases. Since 
the increase in implies more resistance (interference) to the flow, 
it shows that the load capacity and friction force increase as the 
interference due to asperities increases. It can also be noted from 
equations ( 5 . 30 ) and (5-33) that the load capacity and friction force 
increase as or increase. 


If there is no interference of the interacting net-work zone 
in the hydrodynamic region 0 <_ x while in the region Jl. ^ ^ 

(where h = o) 5 the condition of mixed lubrication exist as shown in 
Fig. (5.3), the load capacity and friction force can be calculated 
by following the above procedure and we can obtain the following 
expressions 


U ( £ - £. ) (h - H G ) 

D - \) {h - H G^f 


p £. n ( £ - £. ) 2!I 


( 5 . 34 ) 


(5.35) 


h 


where 




•Ty 


since (h - H a) is an increasing function of M and P. 

V -A. \J 

decreases as M increases, it can be seen from equations (5 #34) and 

X 

( 5 . 35 ) that the load capacity and friction force increase with 
as in the previous case vhich is due to hi^er resistance to the 
fluid flow throu^ interacting net-work zone in region II, £ x £ •<!'). 

It may also seen from equations (5.34) and (5.35) that the load 
qapacily and friction force increase as or increases 

5 . 4 HYDROSIAIIC step bearing UHDER mixed nJBPJCATIQg 

Consider the case of hydrostatic step bearing as shown in 
Pig. ( 5 . 4 ) where the upper plates forms an interacting net work zone 
due to asperities through which the lubricant with viscosity flows. 
In this case, Reynolds equation can be obtained from equation (5-21 ) 
as follows : 
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Integrating equation (5-36) an^ using the ( 
(the flux is constant in this case), we hav 


definition of flow flux 


Q = - 2 tr r P, 


1 dr 


(5.38) 


= - 2ir r P 


2 dr 


Integrating equations 


(5.38) with respect to r and using the boundary 


conditions 


= Pj_ at r - 

Pi = Pq = Ps " ""s 


Pg = 0 at r = rg 


the expressions 


for pressure oau be obtained as follows : 


0 n 

H — j,n — 


Pi " " 27rP^ r^ 


(5.39) 


P2 ~ 2'irP2 


Ps “ ^i ■ 2rP, 


Q S 

21TP2 "" 


^ using 


equations (5.38) and (5.39) the erpression 


for flux can 


b e written as , 


2ir Pj_ 


1_ S, ^ la + 1- iln — 
f: r, F, 


(5.40) 


• u. * +Vi'i Q f'ase can be evaluated by 
The load capacrtjr in this case caii 





W = tr + / ^ 2^ r dr + J ® 2it r P2 dr 


(5.41) 


i s 

which on using equation (5.39) gives, 


Q - 

w = f [ — 


4) (4 - rp 

1 +>- 4 ; ±.] 


1 




(5.42) 


Since functions Pi and P„ are similar to the function P and P as 

I -X! 0 

defined in the equations (5.22) and equation (5.26), and are also 
decreasing function of M, it can be seen from equation (5.42) that for 
a constant flow flux the load capacity increases as M increases i.e. 
as the interference to the flow increases. 


It may also be noted from equation (5.40) that for a given inlet 
pressure the flux decreases as the interference parameter increases. 

Further, as P.^ and P^ also decreases as n.j or increases, it 
is concluded that the load capacity increases as or hg increases for 
given flux but the flow flux decreases for given inlet pressure. 

5.5 SQUEEZE PIIM BEABIgS 

Consider the case of squeeze film through three layers which 
may be applicable to synovial joint functioning under the condition 
of standing, jumping etc. as shown in Pig. (5.5). Under hydrodynamic 
condition of the thickness of the ^novial fluid film is non zero 
but when the joint function under the condition of mixed lubrication, 
h =:0. Ihis condition arises at the fag end of jumping or when falling 








acoidently and creating impact loading condition. In the last 
condition, the cartilage may fractured and ruptured and thus 
affecting the interference parameter In fact, due to the fracture, 

the resistance to flow throir^ the interacting net-work zones would 
decrease and M would decrease. 

(Do see the effect of M , on the characteristics of joint, very 
simply clarly it can be assmed that the lower surface of the 
cartilage is plane while the upper surface may be plane or curved, 
further, it is assumed that under impact loading or near impact 
loading condition, since the cartilages are elastic, the contact zone 
has aconstant film- thickness. 


Keep the above in view and assuming that the thickness of middle 
layer is small but not zero, the equation (5 .20) reduces to the following 
one-dimensional form, as applicable to the present situation, (H^ — 0, 


= H, D, = Uj = 0) 
dp 


Ff -^1 = - V 

dx '■ x2 dx^ 


where 


.0 


M - tanh M 
x2 I2n^ 2n^ ^ jP 


tanh 2M aM + tanh M „ 

+ ( — 2i — ^ ^ 


M 


X 


(5.43) 


(5.44) 


Integrating equation (5.43) and using boundary conditions 


dp 

dx 


= 0 


at X 


2 , 

2 


(5.45) 



(5.46) 


we have, 

dp 

dx 


P 


(x - 4) 


x2 


By considering p = 0 at x = 0 and x = ii , the load capacity of the 
joint can be calculated as follows : 


r <ip / x(x 

W = / (_x— ) dx = 7 / ^ 

0 0 V 


dx 


(5.47) 


The time of squeezing T for a constant load W can be obtained as follows, 
by writing 

dh 




T= -| 1(h) dh 


where H is assmed constant and 


(5.48) 


x(x-4) 
1(h) = / 2 


P 


dx 


(5.49) 


0 ^x2 

Prom the expression of P^ in equation (5.44) and the appendix, it can 
be seen that P^^ decreases as M^, or increase. -Thus in this case 
also it may be noted that the load capacity increases due to decrease 
in P^. Hence, the load capacity and the time of squeezing decrease 
as M , n. , decrease. The decrease in If may be caused due to 
fractured cartilage while decrease in or pg may be caused because 
of hi^ shear rate under impact loading conditions, illso, in the above 
piocess a situation can arise where p^ , the viscosity of the synovial 
fluid film of thickness h may impact be increased because of the presence 



of hyaluronic acid molecules due to boosted lubrication, Dowson, et. 
al. [ 1970 ] . This increase in might compensate its decrease due 
to high shear rate under impact loading condition, 

Bius, it may be concluded that the fractured joint functions 
less effectively than the normal joint. But this mi^t be compensated 
because of the novel structure of the joint provided by nature even 
under severe conditions. 

Under the condition of severe impact loading, the two interacting 
surface net-work zones of the cartilage would be in contact (h = O) 
over the whole region and the two surfaces are almost parallel in this 
region [see Pig. (5.6) ]. In this case, function reduces to the 
following form. 


25:2 h=0 ~ 2n^ 




M - tanh M 




tanh 2x1 aM + tarhM _ 

+ ( £) ( 2 S)2 1 

V 2J.f ^ ^ M ^ 


(5.50) 


X X 

As before, the expressions for load capacity, in this case can be 
obtained as follows : 

!. . 

,3 


/ x:(x - •^) 
W= 7 I = 


?£' 


0 


xo 


12 P 


(5.51) 


xo 


Since the cariiilage are elastic, the thickness of interacting net- 
work zones can decrease because of loading and sQueezing can also 
take place even in this case. Then by considering 
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Y =z ^ "dt ~ ~ ^ f "time of squeezing for a constant load W 

can be obtained as follows : 


2H. 

T = “ / ^ 


6ff 


2H^ 


P 


(5.52) 


xo 


Again, since P increases as M or n. decrease, it-can be seen 
from, the equations (5.51 ) and (5*52) that the load capacity and time 
of squeezing decrease as M or n. decrease due to fractured cartilage 
under the condition of impact loading. 


Thus, a person with a fractured cartilage may not able to stand 
and support his body load due to continuous contact of the cartilage 
surface and quick squeezing of the fluid in the interference zones. 

5.6 COFCIUSIOIE 

A new theoiy of mixed lubrication is presented by considering 
a modified forms of equations of motion in interacting net-work zones 
formed by surface asperities and a generalized form of Reynolds 
equation is derived. 

As particular cases, step slider bearing, hydrostatic bearing 
have been studied. It is shown that the load capacity and friction 
foroe increase as the interference to the flow increases or as 
viscosities of layers increase. 

In the case of ^novial joint, it has been pointed out that the 
load capacity and time of squeezing decrease due to fractured joint 
under impact loading condition. However, this may be compensated by 
boosting effects. 
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APPEEDIX 


- tanh M 

The derivative of the fUnotion — = ^ is given as follows 




m. 


M - tanh M ¥(1 ) 

■X 


X 


where 


.^) = 5M^ - tanh - 3 tanh 


( 5 . 53 ) 


( 5 . 54 ) 


Since, ’i'(O) s 0 and 
d 


dM^ 


'HiHj,) = 2 tanh M^(tanh sech M^) (5.55) 


is positive for all values of >0* Hence ¥(M^) is an increasing 

function of and positive for all values of > 0, Then it ©an be 

M - tanh M 

cohcluded; from equation (5*53)» that — = is a decreasing 




function of for all values of > 0. 


tanh M 


In a similar manner, it can also be shown that 
a decreasing function of M for all values of M > 0, 

X X 


M 


~ is also 


Keeping these in mind, it can be noted from equations (5«22-a), 

(5.26), (5. 37 -a), (5.37-b), (5-44) and (5*50) that the functions P^, P^, 

Pj . P/^. P and P decrease as M increases* It can also be noted that 
1 ^ 2 ' x2 xo 2: 

these functions also decrease as or Pg increases. 

How consider the function, 


tanh M, 


G = [a + 

X 


^ [> 


M ■■ 2M^ 

Hg tanh 21 


(5.56) 


and 





tanh M 

X 

M 


X 


[5,51) 


As before it can be seen from the above equations (5.56) 

(5.57) that and G^ are decreasing functions of M for all values of 

X O w ^ 

M > 0. 


lUrther as M tends to zero. G- and G tend to (a+t)— — ; — — ^ 

X ’ X 0 ^^2 “^1 

and one respectively, igain as tends to “ , G^ or G^ tends to zero, 
Also, 


( a+ i) (ct+l) 


no + an, n, 

‘ a~ + 1 

’'a 


< 1 for > ng 


= 1 for n,^ = ng 


(5.58) 


Since G 3,ncd G‘ are decreasing functions of M , (G. - 6 ) is 

X O X O X 

positive for n, > , and decreases as M increases for all values 

12 X 


of > 0. 


Keeping these in view and noting that, h > H, in the case of 
slider bearing [see Kigs. (5.2) and (5.3)1 and (G^ - G^) <1, it is 
concluded that the functions { h — h( G^ — ^ increase 

as M increases, 

X 
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F 

h 




h,,H2 


3 

A 


h 


PiP^jPg 


Pc 




T 

u,v,w 


“I’M 

ni,Vi 


LISJ OP lOiTAJIOffS 

- friction force 

- filiii-thickness 

- thickness of the interference zones 

- distance between the surfaces and xy plane 

- denoted for ith interface 

- denoted for jth interface 

- bearing length 

- step position 

- fluid film pressure 

- inlet pressure 

- step pressure 

- flux 

- radial coordinate 

- inlet radious 

- outlet radious 

- step position 

- time of squeezing 

- velocities of fluid in x,y and z directions 

- velocities of fluid in z and y directions in the region I 

- velocities of fluid in x and y directions in the region II 

- velocities of fluid in x and y directions in the region III 

- velocities of the surface at z = along x aM y directions. 
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Ui,v. 

u.,v. 

J* 3 
Y 

W 


x,y,z 


velocities of the surface at z = Hg along z and y 
directions. 

velocities at 2-.= H^+ along x and y directions 
velocities at z = along x and y directions 

squeeze velocity 
load capacity 

rectangular coordinate system 

viscosities of lubricant in the regions I, II and III 
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